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ABSTRACT 

An  investigation  is  made  of  the  sound  generated  by  the  impingement  of  a  ventilating 
jet  on  the  gas-water  interface  of  a  ventilated  supercavity.  A  ventilated  supercavity  is  a 
gaseous  envelope  generated  around  an  underwater  vehicle  that  allows  for  order  of  magnitude 
increases  in  vehicle  speeds.  However,  the  hydrodynamic  noise  generated  by  the  supercavity 
can  interfere  with  successful  deployment  of  the  vehicle.  One  of  the  principal  mechanisms  of 
noise  generation  is  believed  to  be  the  impingement  of  the  cavity  ventilating  gas  jets  on  the 
gas-water  cavity  wall.  An  understanding  of  the  acoustic  field  generated  by  this  interaction 
has  been  developed  by  analysis  of  a  series  of  model  problems  which  approximate  the 
geometry  and  physical  mechanisms  involved  in  the  jet-cavity  interaction.  The  first  problem 
is  that  of  a  spherical,  gas-filled  cavity  in  water  whose  surface  is  excited  by  a  planar  ring  of 
axially  projecting  jets  originating  from  the  center  of  the  sphere.  The  second  involves  a  jet 
of  infinitesimal  cross-section  impinging  at  normal  incidence  on  the  gas- water  interface.  The 
final  problem  makes  use  of  a  creeping  mode  diffraction  theory  to  estimate  the  ‘self-noise’ 
produced  by  the  impinging  jets  on  the  solid  nose  of  the  vehicle. 

These  problems  provided  insight  into  the  general  characteristics  of  the  held  and  its 
dependence  on  typical  how  conditions.  To  apply  these  solutions  to  the  problem  of  the 
supercavity,  experimental  studies  were  made  in  collaboration  with  the  Applied  Research 
Laboratory  (ARL)  at  Penn  State  University.  These  included  measurement  of  the  unsteady 
force  exerted  by  a  jet  on  an  interface,  the  results  of  which  have  been  used  to  make  theoretical 
predictions  of  the  self-noise  generated  by  a  model-scale  supercavity  currently  in  use  at  ARL. 
The  final  results  of  this  project  will  aid  current  and  future  naval  research  into  supercavity 
noise  and  methods  for  its  reduction. 
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Chapter  1 


Introduction 


A  ventilated  supercavity  is  a  gaseous  envelope  generated  around  an  underwater  vehi¬ 
cle  that  allows  for  significant  increases  in  vehicle  speed.  However,  the  hydrodynamic 
‘self-noise’  generated  by  the  supercavity  can  interfere  with  successful  deployment  of  the 
vehicle.  In  order  to  reduce  or  eliminate  this  problem  it  is  necessary  to  understand  and 
quantify  the  noise  production  mechanisms.  One  of  the  principal  sources  of  noise  is  be¬ 
lieved  to  be  the  impingement  on  the  supercavity  wall  of  the  cavity-ventilating  gas  jets, 
and  the  outcome  of  this  project  will  be  an  understanding  of  this  mechanism  and  its 
probable  impact  on  the  self-noise  problem. 

1.1  Motivation 

The  successful  development  of  high-speed  supercavitating  underwater  vehicles  is  an  im¬ 
portant  U.S.  Navy  objective.  All  aspects  of  the  technology  are  being  vigorously  promoted 
by  the  Office  of  Naval  Research  and  at  various  naval  laboratories.  The  primary  goal  of 
this  project  is  to  provide  naval  researchers  with  an  understanding  of  the  characteristics 
of  the  sound  generated  by  the  ventilating  jets  of  a  supercavitating  vehicle.  This  will  fa¬ 
cilitate  current  investigations  of  the  overall  sound  generated  by  such  vehicles,  and  permit 
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the  development  of  methods  for  its  control  and  reduction. 

1.2  Background 

1.2.1  Cavitation 

Cavitation  is  the  development  of  a  pocket  of  gas  or  vapor  inside  a  liquid  volume.  It 
is  driven  by  a  drop  in  pressure  to  a  level  below  the  liquid’s  vapor  pressure  causing  it 
to  change  phase  (the  temperature-driven  equivalent  to  this  is  boiling).  Cavitation  can 
occur  in  liquids  either  in  motion  or  at  rest,  and  both  in  the  body  of  the  liquid  or  on 
solid  boundaries.  The  form  of  hydrodynamic  cavitation,  cavitation  in  flowing  liquids, 
of  particular  relevance  to  underwater  vehicles  is  fixed  cavitation,  or  cavitation  which 
develops  attached  to  a  solid  boundary.  The  boundary  geometry  causes  the  detachment 
of  the  liquid  boTindary  layer  from  the  solid  wall,  otherwise  known  as  flow  separation. 
High  velocity  (and  ergo  low  pressure)  flow  allows  for  cavitation  inception  in  the  separation 
region.  Here  small  cavities  coalesce  into  a  large  cavity  attached  to  the  boundary.  Without 
separation  and  this  region  of  relative  calm,  any  small  cavities  which  may  form  are  unable 
to  attach  to  the  boundary  as  they  are  rapidly  convected  away  by  the  liquid  flow  [4,  14,  26]. 

If  favorably  low  pressures  persist,  the  cavity  is  maintained  or  may  continue  to  expand 
downstream  of  the  inception  cite.  Cavities  which  initiate  and  close  on  the  same  solid 
boundary  are  known  as  partial  cavities.  Cavities  which  terminate  downstream  of  the 
boundary  or  body  on  which  they  initiate  are  referred  to  as  supercavities.  The  artificial 
admission  of  gas,  either  into  a  region  of  separated  flow  to  initiate  a  cavity  or  into  a  partial 
cavity  or  supercavity  to  expand  and  maintain  it,  is  called  ventilation.  Ventilated  cavities 
behave  comparably  with  natural  cavities  except  in  the  aft  region  of  gas  exhaustion  and 
fluid  entrainment  and  the  potential  deformation  of  the  cavity  due  to  gravity  (relevant 
only  at  low  flow  velocities)  [4,  14,  26].  The  method  of  gas  exhaustion  as  it  applies  to 
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ventilated  supercavitating  vehicles  is  discussed  below. 


1.2.2  Supercavitating  vehicles 

An  underwater  vehicle  can  attain  high  forward  speeds  by  the  development  of  a  gaseous 
ventilated  supercavity  which  eliminates  skin  friction  and  reduces  drag  by  up  to  90% 
[7,  27].  The  cavity  is  initiated  downstream  of  a  specially  designed  ‘cavitator’  on  the 
vehicle  nose,  ft  is  expanded  to  contain  the  entire  aft  section  of  the  vehicle,  except  for 
control  surfaces  used  for  guidance,  by  the  injection  of  gas  from  sources  just  aft  of  the 
cavitator  (Fig.  1.1).  The  cavity  can  be  maintained  in  a  stable  configuration  by  carefully 
controlling  the  gas  injection  rate  to  avoid  overpressures  and  instabilities  that  can  cause 
the  cavity  to  pulsate  [14,  27,  41,  43].  Under  stable  conditions,  involving  high  speed 
motion  through  the  water,  the  gas  exhausts  from  the  rear  end  of  the  cavity  by  the 
quasi-periodic  shedding  of  ring  vorticies  with  gaseous  cores.  At  low  speeds  (i.e.  at  low 
Froude  numbers  that  are  not  necessarily  relevant  in  the  present  context)  the  supercavity 
resembles  a  hot  ’plume’  and  gas  escapes  through  two  trailing,  hollow  vortex  tubes  formed 
by  buoyancy-induced  plume  bifurcation  [14,  39]. 


Mean  Flow  i 


Figure  1.1.  Schematic  experimental  supercavity 


The  supercavity  is  a  significant  source  of  aerodynamic  sound  that  can  interfere  with 
the  underwater  vehicle’s  guidance  systems,  inhibiting  successful  deployment  [6,  18,  19,  32] 
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(see  Appendices  A  and  B  for  a  brief  review  of  concepts  from  acoustic  theory  relevant  to 
this  work).  Turbulence  in  the  aqueous  boundary  layer  approaching  the  wetted  trailing 
edge  of  the  cavitator  from  the  nose  generates  sound  and  hydrodynamic  pressure  fluctua¬ 
tions  as  it  convects  across  the  edge  [19,  21,  48].  To  analyze  sources  further  downstream 
it  is  reasonable  in  a  first  approximation  to  neglect  pressure  fluctuations  within  the  cav¬ 
ity  because  of  the  vast  difference  in  mass  densities  of  the  gas  and  water,  in  which  case 
the  gas-water  interface  can  be  treated  as  a  pressure-release  surface.  Then  the  sound 
generated  by  turbulence  quadrupoles  in  the  flow  adjacent  to  the  cavity  edge  should  be 
relatively  unimportant,  because  the  pressure-release  interface  causes  the  sound  pressure 
to  vary  as  pwv2M 3,  where  the  Mach  number  M  =  v/cw  -C  1,  v  being  the  flow  velocity 
and  pw,  cw  are  respectively  the  mean  density  and  sound  speed  in  water  [12],  Bubbles 
and  water  droplets  in  the  break-up  region  of  the  cavity  far  downstream  are  respectively 
equivalent  to  monopole  and  dipole  sources;  except  at  extremely  high  frequencies  their 
importance  is  again  greatly  reduced  by  the  proximity  of  the  pressure-release  cavity. 

Observation  suggests  that  the  acoustic  noise  is  dominated  over  a  broad  range  of  fre¬ 
quencies  by  the  unsteady  impingement  of  the  ventilating  jets  on  the  cavity  wall  [48].  The 
jets  ripple  the  gas- water  interface,  producing  unsteady  surface  forces  and  a  distributed 
acoustic  surface  source  of  dipole  type.  In  order  to  investigate  the  noise  from  this  source 
independently  from  other  mechanisms  at  work  in  and  around  the  supercavity,  a  theoret¬ 
ical  description  of  the  impingement  mechanism  as  an  acoustic  source  is  required. 

To  facilitate  the  development  of  these  mathematics,  we  will  consider  a  slightly  modi¬ 
fied  gas  injection  system.  On  typical  ventilated  supercavitating  vehicles,  a  gas  deflector 
is  fitted  over  the  gas  injection  site,  as  in  Figure  1.2.  This  redirects  the  jets  to  impact  the 
cavity  at  shallow  incidence.  We  will  choose  instead  to  consider  that  the  gas  deflector  is 
absent  (as  in  Fig.  1.3),  such  that  the  jets  impact  the  cavity  at  normal  incidence.  This 
configuration  allows  for  much  simpler  geometric  approximations  of  the  jet-cavity  inter- 
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action  to  be  used  in  the  theoretical  work.  Furthermore,  the  removal  of  the  gas  deflectors 
from  the  experimental  model  will  promote  sound  production  via  the  impingement  mech¬ 
anism,  expediting  the  in- water  acoustic  surveys  of  the  supercavity  to  be  performed  by 
naval  researchers. 


Figure  1.2.  Experimental  vehicle  inside  ARL  305  mm  water  tunnel 


mean  flow  jet  impact 


Figure  1.3.  Experimental  vehicle  inside  ARL  305  mm  water  tunnel  with  gas  deflectors 
removed  and  jets  impinging  at  normal  incidence  on  the  cavity  interface 
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1.2.3  Discussion  of  existing  literature 

A  literature  survey  revealed  that  the  dynamics  of  jet  impingement  are  well  explored. 
Extensive  experimental  and  theoretical  studies  of  impingement  on  both  rigid  and  free 
surfaces  have  been  made  (see  [17,  36,  46]  and  their  bibliographies  as  examples  thereof). 
A  significant  amount  of  attention  has  been  payed  to  the  heat  transfer  problem  of  jet 
impingement  on  a  plate  due  to  its  prevalence  in  manufacturing  processes,  however  such 
work  is  of  no  relevance  for  the  supercavity  problem.  The  acoustics  of  a  jet  impinging  on  a 
plate  have  likewise  been  explored  [28,  34,  40],  though  almost  exclusively  experimentally. 
Furthermore,  these  studies  focus  on  mechanisms  of  noise  generation  (such  as  shed  vor- 
ticity  and  turbulence)  other  than  the  unsteady  motion  of  a  free  surface.  No  theoretical 
description  of  the  sound  generated  by  jet  impingement  on  a  gas-water  interface  exists; 
it  is  necessary  to  develop  this  formulation  for  the  acoustic  field  using  the  equations  and 
methods  fundamental  to  the  study  of  hydroacoustics.  The  few  existing  studies  which 
may  then  prove  relevant  to  the  supercavity  problem  are  those  which  provide  information 
about  the  unsteady  pressure  exerted  by  the  jet  on  an  impingement  surface  (as  in  Strong 
et  al.  [45]).  These  experimental  results  may  serve  as  a  ‘check’  of  empirical  studies  of  the 
supercavity  jets  and  of  the  generality  of  the  mathematical  formulation  of  the  generated 
sound. 


1.3  Approach 

The  problem  was  investigated  by  consideration  of  a  series  of  model  problems  and  related 
experiments  involving  predictions  made  on  the  basis  of  approximations  to  the  geometry 
and  physical  mechanisms  governing  the  jet-cavity  interaction.  The  simplest  approxima¬ 
tion  involves  a  ventilated  supercavity  in  the  form  of  a  spherical,  gas-filled  cavity  in  water 
whose  interior  surface  is  excited  by  a  planar  ring  of  axially  projecting  jets  originating 
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from  the  center  of  the  sphere.  This  problem  was  not  designed  to  predict  the  acoustic 
field  generated  by  the  supercavity  jets;  it  was  a  convenient  starting  point  which  provided 
general  insight  into  the  behavior  of  the  held  and  the  source  mechanisms  at  work  in  the 
supercavity.  The  solution  of  the  sphere  problem  yielded  insight  into  the  directivity  of 
the  sound  radiating  into  the  water  and  the  likely  contributions  from  cavity  resonances; 
it  also  enabled  the  exploration  of  the  importance  of  several  damping  mechanisms. 

The  second  problem  examined  was  the  production  of  sound  by  a  jet  of  infinitesimal 
cross-section  impinging  at  normal  incidence  on  a  planar  gas-water  interface.  Its  solution 
was  formulated  in  terms  of  the  force  exerted  by  the  jet  on  the  cavity  wall.  The  resulting 
prediction  of  the  acoustic  held  helps  in  understanding  the  more  complex  issues  involved 
in  the  related  problem  of  jet  impingement  on  a  cylindrical  cavity. 

To  apply  these  general  solutions  to  the  problem  of  the  supercavity  an  empirical  knowl¬ 
edge  of  the  gas  injection  jets  was  required.  This  knowledge  was  gained  through  experi¬ 
ments  made  in  collaboration  with  the  Applied  Research  Laboratory  (ARL)  at  Penn  State 
University.  These  experiments  included  measurement  of  the  unsteady  force  exerted  by  a 
jet  on  an  interface,  the  results  of  which  were  then  used  to  make  theoretical  predictions 
of  the  self-noise  generated  by  a  model-scale  supercavity  in  use  at  ARL. 

The  results  were  used  finally  to  estimate  the  contribution  to  the  self-noise  at  the 
vehicle  nose  from  the  jet  impingement  source.  This  was  accomplished  by  use  of  a  creeping 
mode  diffraction  theory  to  model  propagation  of  sound  over  the  rounded  nose  region.  To 
do  this  the  supercavity  was  modeled  as  a  semi-infinite  cylindrical  cavity  to  the  rear  of  a 
solid,  elliptical  nose  section  with  a  prescribed  surface  impedance.  The  cavity  is  excited 
by  gas-jet  impingement  downstream  of  the  cavitator  which  generates  a  separately  derived 
pressure  field.  This  pressure  held  was  used  in  the  diffraction  theory  to  predict  the  sound 
at  the  nose  and  thereby  permit  more  definitive  conclusions  to  be  drawn  regarding  the 
presence  of  noise  in  the  vicinity  of  the  guidance  system  at  the  nose. 
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1.4  Outline 

Chapter  2  describes  the  general  problem  of  a  spherical  ventilated  cavity.  The  sound 
produced  by  jet  impingement  on  a  planar  gas-water  interface  is  explored  mathematically 
in  Chapter  3.  The  experiments  measuring  the  force  exerted  by  the  scale-model  vehi¬ 
cle’s  gas  injection  jets  are  detailed  in  Chapter  4,  as  are  preliminary  predictions  for  the 
impingement  generated  sound.  The  impingement  noise  heard  at  the  nose  is  determined 
in  Chapter  5.  Chapter  6  summarizes  the  body  of  work  and  its  findings  and  presents 
recommendations  for  future  research  into  ventilated  supercavity  acoustics. 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


Chapter  2 


Sound  generated  by  a  jet-excited 
spherical  cavity 


This  chapter  examines  the  general  behavior  of  the  acoustic  held  generated  by  the  im¬ 
pinging  jets.  A  spherical  cavity  in  which  gas  enters  from  the  center  in  a  ring  of  radially 
projecting  jets  is  considered.  This  highly  simplified  geometry  allows  for  general  deduc¬ 
tions  about  the  generated  sound  without  the  difficulty  implicit  in  a  more  complicated, 
realistic  cavity  geometry.  The  mathematical  problem  is  formally  solved  by  series  expan¬ 
sion  in  Section  2.1.  Section  2.2  discusses  low  frequency  cavity  resonances.  High  frequency 
sound  is  addressed  in  Section  2.3.  Specific  results  and  conclusions  are  drawn  by  choosing 
the  volume  of  the  cavity  to  be  approximately  the  same  as  that  of  a  typical  experimental 
supercavity.  The  work  presented  in  this  chapter  has  also  appeared  in  Journal  of  Sound 
and  Vibration  [13]. 
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2.1  The  acoustic  problem 

2.1.1  Problem  formulation 


Figure  2.1.  Configuration  of  the  spherical  cavity 

Consider  a  spherical  cavity  of  undisturbed  radius  a  in  water  of  mean  density  pw  and 
sound  speed  cw  which  is  at  rest  at  infinity.  Let  the  origin  O  be  at  the  center  of  the  cavity 
(Fig. 2.1)  and  introduce  spherical  polar  coordinates  (r,  6,  (ft)  (0  <  6  <  n,  0  <  </>  <  27t), 
where  the  latitude  6  is  measured  from  the  positive  x  direction  in  the  figure.  Gas  enters 
the  cavity  via  a  thin,  axisymmetric  planar  jet  from  a  source  at  O,  and  impinges  normally 
on  the  cavity  interface  S  along  the  great  circle  9  —  |.  It  is  assumed  that  the  mean 
volume  of  the  cavity  and  the  mean  density  pQ  of  the  gas  in  the  cavity  are  maintained 
constant  by  the  steady  exhausting  of  gas  into  the  water  from  the  ‘rear  end’  of  the  cavity, 
say  the  vicinity  of  the  point  A  in  Fig.  2.1. 

Turbulence  fluctuations  in  the  jet  are  assumed  to  be  axisymmetric  and  to  exert  an 
unsteady  radial  force  on  the  cavity  interface  equal  to  F(t)  per  unit  length  of  the  circle 
of  impingement  on  S.  The  definition  of  F(t)  is  made  precise  by  requiring  it  to  equal  the 
force  that  would  be  exerted  on  the  interface  by  the  jet  when  the  interface  is  assumed  to 
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be  rigid.  The  forced  motion  of  the  interface  produces  pressure  fluctuations  p(r,9,t )  that 
satisfy 


\  c2  dt2 


(1_  (P_ 
V4  dt 2 


0,  r  <  a, 

0,  r  >  a, 


(2.1) 


respectively  within  the  cavity  and  in  the  water,  where  ca  is  the  mean  speed  of  sound  in 
the  gas. 

The  analytical  problem  is  simplified  by  introducing  the  Fourier  decomposition  F(t)  = 
F0(uj)e^lujt  dui  and  considering  first  the  solution  for  the  case  where  the  cavity  is 
excited  by  a  time-harmonic  radial  force  F0(u)  per  unit  length.  The  solution  of  the  time- 
dependent  problem  can  subsequently  be  found  by  application  of  the  integral  operator 
fZ(  ■  )e~tujt  duj.  Then  time- harmonic  pressure  p(r,  9,  u)  satisfies 


(V2  +  k2)  p  =  0,  where 


k  =  ka  =  u/co,  r  <  a 
k  =  kw  =  uj / cw,  r  >  a 


(2.2) 


The  solution  in  the  outer  region  r  >  a  must  exhibit  outgoing  wave  behavior.  The  two 
solutions  in  r  >  a  are  related  by  the  conditions  that  the  pressure  and  normal  displacement 
are  continuous  at  the  interface.  Because  the  force  per  unit  area  applied  to  the  interface 
by  the  jet  is  ( F0/a)5  (9  —  §),  the  pressure  will  be  continuous  at  r  =  a  provided 


p(a  +  0,  9,  to)  =  p(a 


0, 9 ,  oS)  + 


f )  ,  0  <  9  <  7T. 


(2.3) 


Similarly,  because  F0  corresponds  to  the  force  exerted  on  a  rigid  interface,  continuity  of 


11 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


normal  displacement  is  satisfied  by  requiring 


-(?r) 

Pw  \drjr=a+ o 


1(tt)  • 

Po  \drjr=a_ o 


0  <  9  <  7T, 


(2.4) 


because  velocity  fluctuations  within  the  jet  at  the  interface  are  balanced  by  the  unsteady 
pressure  gradient  in  the  jet. 


2.1.2  Solution  in  the  general  case 

Within  the  cavity  the  time-harmonic  solution  that  remains  finite  at  r  =  0  has  the 
expansion  [18,  44] 

OO 

P  =  yAnjn(k0r)Pn(cos9),  r  <  a,  (2.5) 

n= 0 

where  jn  is  a  spherical  Bessel  function  of  the  first  kind,  Pn  denotes  the  Legendre  poly¬ 
nomial  of  order  n  [1],  and  the  coefficients  An  are  to  be  determined.  In  the  water  the 
pressure  is  expanded  in  terms  of  outgoing  radiating  waves  in  the  form 

OO 

P  =  Bnh^)(kwr)Pn(cos9),  r  >  a,  (2.6) 

n= 0 

where  the  Bn  are  constants  and  is  the  spherical  Hankel  function  of  the  first  kind  [1], 
The  coefficients  An,  Bn  are  calculated  from  conditions  (2.3)  and  (2.4).  The  first 
supplies 

OO 

(Bnh^\kwa)  -  Anjn(k0a))  Pn(cos  9)  =  —5  (0  -  f)  ,  0  <  9  <  vr,  (2.7) 

'  'CL 

n= 0 

which,  by  means  of  the  orthogonality  relation  fQ  Pn( cos  9)Pm( cos  9)  sin  9  d9  =  2Snm/ (2n+ 
1)  and  the  formula  P„(0)  =  cos(I|?:)r  (f  +  |)  j 'JkT  +  l)  [1],  reduces  to 
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Bnhn  ( kwa )  Anjn  (fc0Cl) 


F0  (2n  +  1) 
a  2^/7? 


cos 


r(t  +  |) 

r(f  +  i)' 


(2.8) 


Condition  (2.4)  yields 


B"k"'K'\Ka)  -  ^'(ha)  =  0. 


Pu 


Po 


(2.9) 


where  the  prime  denotes  differentiation  with  respect  to  the  argument.  Hence 

(2 n  +  1)  cog  (nn\  ^  (2  +  2)  _ _ 

2^/tt  2  r  (I  +  l)  ajn(k0a)(Zn  -  1) 

(2n  + 1)  r(|  +  |)  F0Zn 

2^  '2T(2  +  l)  afc!,,|(i„,a)(Z„-l) 

where 

^  =  pwcw  jn(k0a)  h^(kwa) 

Po^o  3ni)^o&)  fan  (A/i yd) 

Now  cos(^f)  =  0  when  n  is  odd.  Therefore,  by  replacing  n  in  the  above  expressions 

by  2 n,  using  the  result  [1] 


(2.10) 


(2.11) 


A„  — 


= 


r (n  +  1)  _  rl  •  3  ■  5...(2ra  -  1)  _  r  (2n)! 
T(n  +  1)  “  V?r  2^!  =  V7T(2 -n!)2’ 


and  substituting  into  the  expansions  (2.5)  and  (2.6),  we  deduce  the  following  desired 
representations  of  the  sound  within  the  cavity  and  radiated  into  the  water: 


/?  00 

pM,w)  = 

n=0 

(4n  +  l)(2n)! 

!  P2n( c 

(2nn!)2 

(^2n 

00 

=  2a 

n= 0 

(4n  +  l)(2n)! 

(2nn!)2 

(3, 

r  <  a, 


(Z2n-l)h^(kwa)  ’ 


r  >  a. 


(2.12) 


(2.13) 
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2.1.3  Special  limiting  cases 


It  is  also  useful  to  note  the  corresponding  forms  of  equations  (2.12),  (2.13)  when  (i)  the 
water  is  regarded  as  incompressible  (i.e.  when  cw  — >  oo),  and  (ii)  there  is  no  gas  in  the 
cavity  (pa  — >■  0). 

In  case  (i)  there  can  be  no  damping  of  the  sound  in  the  cavity  by  radiation  losses, 
and  we  find 


p{r,9,u) 


Fo  vy  (4n  +  l)(2n)!  P2n{cos0)j2n(kor) 
2a  ^  (2nrd)2  {Z'2n  -  l)j2n{k0a)  ’ 


(2.14) 


—  Vi  p"(4n+1)(2n)!  Z2nP2n(cOsg)  /a\2n+l 
2a  J  (2nn!)2  (Z'n-1)  \rJ  ’ 


(2.15) 


where 

f^l  _  PwkoQ'  J 

”  Po(n  +  1)  jn(k0a) ' 


(2.16) 


In  case  (ii),  where  the  cavity  is  treated  as  a  vacuum,  Zn  — >  oo  and  there  are  no 
pressure  fluctuations  within  the  cavity  (p  =  0  for  r  <  a).  Then  expression  (2.13)  reduces 
to 


P(r,9,w) 


Fo  yy  ,n(4u  +  l)(2ra)!  P2n(cosd)h^(kwr) 

2aho  h*nM 


(2.17) 


2.1.4  Suppression  of  coherent  interior  modes 

Just  inside  the  cavity  interface  the  Bessel  function  jn(k0r )  =  \  (hn\k0r)  +  hn\k0r )^  in 
the  expansion  (2.5)  of  the  cavity  pressure  can  be  interpreted  as  representing  a  compo¬ 
nent  traveling  wave  ~  hn\k0r )  incident  on  the  interface  from  within  the  cavity,  and  a 
component  ~  hn\kar )  reflected  from  the  interface.  Waves  reflected  from  the  interface 
will  subsequently  impinge  again  on  the  interface,  possibly  leading  to  the  establishment 
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of  coherent  resonant  oscillations.  At  high  frequencies  scattering  by  small  scale  irregulari¬ 
ties  on  the  interface  will  induce  incoherence  in  the  interior  reflected  wave  field  tending  to 
suppress  the  growth  of  coherent  resonances.  By  discarding  the  incident  coherent  compo¬ 
nent  hn\k0r )  of  the  Bessel  function  we  can  therefore  hope  to  mimic  the  effect  of  random 
interface  scattering. 

This  is  equivalent  to  repeating  the  analysis  leading  to  the  predictions  (2.12)  and 
(2.13)  with  jn{k0r)  replaced  throughout  by  hn\kQr ),  and  yields 


P(r,0,w)  = 


,(4n  +  l)(2n)!  P2n (cos 9)h(^{k0r) 


n=0 


(2nn!)2  (Z2n 


1  )h?n(k0a) 


r  <  a, 


(2.18) 


Fo  \- .  ,„(4ti+ 1)(2ti)1  Z2nP2n(c os 

2a2J-  (2"n!)2  (Z2„  -  l)h%(k„a) 


where 

^  _  pwcwh^'(k0a)  h^(kwa) 

n  PoCo  h$(k0a)  hn)f(kwa)  1  J 

2.2  Low  frequency  resonances 

The  spectrum  of  sound  radiated  to  large  distances  from  a  cavity  tends  to  be  dominated  by 
low  frequencies,  comparable  to  the  fundamental  Minnaert  volumetric  pulsation  frequency 
[6,  7].  This  and  other  resonant  oscillations  of  the  cavity  are  excited  by  the  jet.  In  this 
section  we  examine  this  excitation.  Predictions  will  be  illustrated  for  the  case  of  an 
air-filled  cavity  in  water  at  standard  temperature  and  pressure  (STP)  when  the  cavity 
has  radius  a  =  7.0  cm,  which  corresponds  to  a  cavity  volume  approximately  the  same  as 
that  of  a  typical  small  scale  experimental  supercavity. 
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2.2.1  Monopole  resonances 

The  monopole  is  nominally  the  most  efficient  source,  but  it  is  relevant  only  at  very  low 
frequencies.  Volumetric  pulsations  are  then  analogous  to  oscillations  of  a  mass-spring 
system  in  which  the  moving  ‘mass’  is  the  water  displaced  radially  by  the  cavity,  and  the 
‘spring’  stiffness  is  furnished  by  the  compressibility  of  the  contained  gas.  The  compress¬ 
ibility  of  the  water  makes  a  negligible  contribution  to  the  gross  motion  of  the  cavity  at 
low  frequencies,  when  the  acoustic  wavelength  is  very  much  larger  than  the  radius  a  of 
the  cavity.  The  resonance  frequencies  can  therefore  be  calculated  by  ignoring  the  aqueous 
compressibility,  which  is  responsible  only  for  relatively  weak  radiation-damping  of  the 
oscillations.  This  is  the  only  damping  mechanism  available  to  an  ideal  spherical  cavity 
when  viscous  and  thermal  effects  are  ignored  [11,  30].  However,  our  model  implicitly 
assumes  that  the  gas  content  of  the  cavity  is  maintained  roughly  constant  by  the  steady 
exhausting  of  the  gas  into  the  water,  and  we  shall  estimate  below  the  damping  arising 
from  this. 

According  to  equations  (2.14),  (2.15)  resonant  oscillations  of  a  mode  of  order  2 n  are 
determined  by  the  zeros  of  ( Z'2n(k0a )  —  1)  j2n(k0a)  =  0  when  the  water  is  regarded  as 
incompressible,  where  Z'2n  is  defined  as  in  equation  (2.16).  Using  the  relation  j'n(z)  = 
jn-i(z)  —  ^~~jn (2) ,  the  resonance  condition  is  equivalent  to 


Po 

Pw 


(2.21) 


For  the  monopole  n  =  0,  (2.21)  reduces  to  [1] 


(2.22) 


When  pQ/pw  <C  1  the  smallest  positive  root  of  this  equation  is  kQa  ~  \/3p0/pw ,  which  cor¬ 


responds  to  volumetric  oscillations  at  the  Minnaert  radian  frequency  u  =  \J3 p0c20/ pwa2 
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[11,  30]  for  which  the  pressure  does  not  vary  with  position  in  the  cavity. 

For  our  model  cavity  with  a  =  7cm  the  Minnaert  frequency  (ju/2n  ~  47 Hz.  But 
this  mode  is  actually  atypical  of  monopole  resonances,  because  it  corresponds  to  the 
special  case  in  which  the  acoustic  wavelength  is  much  larger  than  the  diameter  of  the 
cavity.  Higher  frequency  volumetric  resonances  (determined  by  the  larger  roots  of  equa¬ 
tion  (2.22))  involve  pressure  fluctuations  within  the  cavity  that  oscillate  between  positive 
and  negative  values  with  distance  from  the  center.  The  first  ten  monopole  modes  are 
listed  in  Table  2.1:  the  wavelengths  27 r/k0  of  the  higher  order  modes  ( N  >  1)  are  always 
smaller  than  the  cavity  diameter. 


mode  number  1V1234  5  6  7  8 

kQa  0.06  4.49  7.73  10.90  14.07  17.22  20.37  23.52 

frequency  Hz  47  3,474  5,972  8,429  10,874  13,312  15,748  18,181 


Table  2.1.  Cavity  monopole  wavenumbers  and  frequencies  at  STP  (from  Eq.  (2.22)). 


The  TVth  order  monopole  resonance  wavenumber  kQ  is  well  approximated  at  higher  fre¬ 
quencies  by  the  formula 


{k0a)N 


(2JV  -  1)7T  r 


"1(1  Pol  Pw ) 
(2 N  -  1  )%2 


for  N  >  4. 


(2.23) 


2.2.2  Influence  of  dissipation 

The  amplitudes  of  resonant  modes  excited  by  the  jet  are  governed  by  losses  produced  by 
acoustic  radiation,  and  by  viscous  and  thermal  diffusion  [7,  11],  and  also  by  hydrodynamic 
convection  of  vibrational  energy  from  the  cavity  by  the  steady  outflow  of  gas  into  the 
water  [18].  These  mechanisms  typically  produce  small  corrections  in  the  values  of  the 
resonance  frequencies  predicted  in  the  absence  of  dissipation.  The  first  order  effect  is  that 
each  undamped  resonant  frequency  u>,  say,  is  replaced  by  u;  ~  cD  —  ie,  where  e  =  e(u>)  >  0 
is  a  small  imaginary  part  that  would  cause  the  corresponding  unforced  resonance  to 
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decay  in  amplitude  like  e_et.  According  to  Devin  [11],  the  separate  contributions  to  e 
from  radiation,  viscous  and  thermal  damping  are  all  of  comparable  magnitudes.  We  can 
evaluate  the  effect  of  radiation  from  the  solution  (2.12),  (2.13)  which  takes  account  of 
the  compressibility  of  the  water,  i.e.  from  the  solutions  of  Z2n(k0a)  —  1  —  0: 


PwCwj2n(k0a)  h(2n(kwa) 
PoCq  j2n(k0&)  h<2n  ( kwQ, ) 


1,  where  k0 


-  k 

Co 


UJ 


(2.24) 


It  will  be  sufficient  to  consider  the  damping  of  the  monopole  modes  (n  =  0),  and 
in  particular  the  dominant  pulsational  mode,  for  which  ua/c0  =  kQa  ~  0.06  (N  =  1  in 
Table  2.1).  By  setting  n  =  0  in  equation  (2.24)  we  derive  the  analogue  of  the  undamped 
frequency  equation  (2.22) 


kQa  — 


—  +  ikQa  °  j  tan (k0a)  =  0, 
Pw  Pw^w J 


(2.25) 


from  which  the  compressible  correction  to  the  Minnaert  approximation  is  readily  deduced 
to  be: 

kQa  ~  kQa  [  1  —  i  °  ]  ,  where  kQa  =  , (2.26) 

Y  2 cw  J  y  pw 

Therefore  the  minimum  complex  frequency  of  volumetric  resonances  is  given  by 


e  = 


c20(k0a) 
2  cwa 


2 


(2.27) 


where  k0a  ~  0.06  (Table  2.1). 

This  estimate  of  the  radiation  damping  coefficient  e  should  be  characteristic  of  the 
damping  also  expected  from  viscous  and  thermal  losses.  However,  it  is  important  that 
it  be  compared  with  the  hydrodynamic  damping  produced  by  the  efflux  of  gas  from  the 
‘rear’  of  the  cavity.  This  is  determined  by  the  following  argument. 

Let  S  denote  the  mean  energy  density  (per  unit  volume)  of  the  resonance  mode 
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within  the  cavity.  The  rate  at  which  energy  of  this  mode  is  removed  from  the  cavity  by 
convection  in  the  gas  exhausting  into  the  water  is  Q£ ,  where  Q  is  the  volume  velocity 
of  the  exhausting  gas,  which  is  the  same  as  the  volume  inflow  rate  in  the  jet.  Hence,  if 
V  =  |7r a3  is  the  cavity  volume,  the  net  monopole  energy  VS  within  the  cavity  decays  at 
a  rate  determined  by 


d£  _  _Q£_ 
~dt  ~  ~~Vr' 


(2.28) 


But  £  ~  \p\2  /  p0c2  and  d\p\2/dt  =  —  2e|p|2,  where  e  is  the  corresponding  imaginary  part 
of  the  resonance  frequency  associated  with  hydrodynamic  damping,  i.e.  e  =  Qj 2V. 
Hence,  from  (2.27),  (2.28) 


^hydrodynamic 

^radiation 


Q/V 


cl(k0a)2  /  cwa 


(2.29) 


The  typical  gas  inflow  rate  Q  ~  0.01  m3/s  for  our  notional  experimental  supercavity 
whose  volume  is  the  same  as  that  of  a  sphere  of  radius  a  =  7  cm,  Taking  k0a  =  0.06  and 
using  STP  values  for  the  other  quantities  in  (2.29),  we  deduce  that 


^hydrodynamic  n 

-  ~  U.U /, 

^radiation 

and  therefore  that  the  steady  efflux  of  gas  into  the  water  has  no  significant  impact  on 
damping. 

This  conclusion  is  doubly  true  for  higher  frequency  modes  which  radiate  more  freely 
into  the  water.  The  implication  is  that  the  acoustic  properties  of  the  cavity  are  not  sig¬ 
nificantly  influenced  by  the  gaseous  exhaust.  This  is  evidently  a  consequence  of  pressure 
continuity  across  the  cavity  interface  and  in  the  exhaust  flow. 
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2.2.3  Quadrupole  resonances 

The  case  n  —  1  in  equation  (2.13)  represents  a  radiation  field  of  quadrupole  directiv¬ 
ity.  The  corresponding  resonance  frequencies  when  the  compressibility  of  the  water  is 
neglected  are  determined  by  setting  n  —  1  in  equation  (2.21),  i.e.  by  the  roots  of 

( k0a )3  —  9 kQa  ( 1  —  —  J  +  ( 9  —  4 (k0a)2  —  ^-^(3  —  (fc0a)2N)  ]  tan(/cDa)  =  0.  (2.30) 

V  PwJ  V  Pw  '  '/ 

The  length  scales  of  these  wave  modes  within  the  cavity  are  no  larger  than  the  cavity 
diameter.  Indeed  it  is  readily  confirmed  by  expansion  of  the  equation  in  powers  kQa  that 
there  are  no  solutions  satisfying  \kQa\  -C  1.  Numerical  solution  of  equation  (2.30)  yields 
the  results  in  Table  2.2,  where  the  frequencies  tabulated  in  the  third  line  are  for  a  7 cm 
spherical  air  cavity. 

mode  number  1V1234  5  6  7  8 

kQa  3.34  7.29  10.61  13.85  17.04  20.22  23.39  26.55 

frequency  Hz  2,585  5,636  8,205  10,704  13,175  15,632  18,082  20,526 

Table  2.2.  Cavity  quadrupole  wavenumbers  and  frequencies  at  STP  (from  Eq.  (2.30)) 

Similar  analyses  of  the  resonances  can  be  performed  for  the  higher  order  modes  {n  >  1 

in  equations  (2.12),  (2.13)),  but  it  is  not  necessary  to  present  details  here.  It  is  clear  that 

the  principal  low  frequency  resonances  occur  for  kQa  >  1,  the  only  exception  being  the 
special  case  of  the  Minnaert  monopole.  Experiments  on  cavitating  jets  all  indicate  that 
the  far  field  acoustic  spectrum  peaks  in  the  neighborhood  of  a  frequency  comparable  to 
the  Minnaert  frequency  [6,  7]  and  typically  decays  like  u;~2  at  higher  frequencies.  Our 
results  indicate  a  similar  behavior  for  the  cavity  (the  corresponding  continuum  spectrum 
being  formed  by  the  aggregate  contributions  from  the  cavity  resonances)  although  the 
precise  details  for  the  spherical  cavity  and  for  a  practical  ventilated  supercavity  depend 
on  the  frequency  dependence  of  the  appropriate  jet  forcing  function  F0(uj). 
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2.3  Directivity  of  the  sound 

The  acoustic  far  field  in  the  water  satisfies  kwr  1.  In  this  limit  the  asymptotic 
approximation  h^(z)  ~  (— i)n+1elz/z  permits  the  reduction  of  equation  (2.13)  to  the 
form 


p(r,  9 ,  to)  ~ 


ijp  e-iu(t-r/cw) 
kwar 


OO 


£ 

n= 0 


An  Z2nP 2n(cC)S  $) 
(Z2n  -  1  )h^n{kway 


oo, 


(2.31) 


where  An  =  (4 n  +  l)(2n)!/{2(2"+1)(n!)2}.  Therefore 


p(r,  9,  u) 

2  1 

\  ^  AnArnZ2nZ2mP 2n(cOS$).F> 27n(cOS$) 

F0/a 

~  ( kwr )2 

Jbo  ~  1)(z2m  -  1  )h^(kwa)h{2Z(kwa) 

(2.32) 


This  formula  determines  the  directivity  (oc  \p(r,  9,uj)\2)  of  the  sound  radiated  into 
the  water.  The  directivity  is  plotted  as  a  function  of  9  in  Fig.  2.2  (on  a  polar  plot,  each 
curve  being  normalized  with  respect  to  its  maximum  value)  for  a  cavity  filled  with  air 
at  standard  temperature  and  pressure  when  kwa  =  0.1,  1,  10,  100.  The  field  shape  is 
spherically  symmetric  at  very  low  frequencies  ( kwa  <  0.1)  characteristic  of  the  Minnaert 
monopole.  As  the  frequency  increases  (for  kwa  greater  than  about  10,  i.e.  for  frequencies 
exceeding  38,600  Hz  for  a  spherical  cavity  of  radius  7  cm)  the  sound  exhibits  a  progressive 
tendency  to  radiate  preferentially  in  the  side-line  direction  with  the  directivity  of  the 
principal  side  lobe  being  approximately  sin2  9,  i.e.  that  of  an  axisymmetric  dipole  with 
axis  along  the  radial  direction  of  the  jet. 

This  dipole  radiation  is  accompanied  by  narrow  beams  radiated  in  the  fore  and  aft 
directions  (9  =  0,  n).  The  presence  of  these  beams  appears  to  be  an  artifact  of  the 
spherical  geometry,  and  would  likely  be  absent  for  the  more  realistic  geometries,  such 
as  that  of  the  ‘cigar’  shaped  cavity  illustrated  in  Figure  1.1,  because  the  interface  is 
effectively  equivalent  to  a  ‘pressure  release’  surface  along  which  sound  cannot  propagate. 
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Figure  2.2.  Directivity  of  the  jet-excited  sound  for  the  general  case  of  an  air-hlled 
spherical  cavity  at  STP  determined  by  Eq.  (2.32)  when  kwa  =  0.1,  1,  10,  100.  Each 
curve  is  normalized  with  respect  to  its  peak  value. 


This  is  illustrated  by  reference  to  the  exact  analytical  solution  for  the  case  of  a  two- 
dimensional  jet  impinging  on  a  plane  interface  (Figure  2.3). 


Figure  2.3.  The  problem  of  the  sound  produced  by  a  two-dimensional  jet  impinging  on 
a  nominally  plane  interface 


For  the  simplest  condition  in  which  the  pressure  in  the  cavity  (the  lower  region  in 
Figure  2.3)  is  negligible,  the  pressure  radiated  into  the  water  consists  entirely  of  the 


22 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


dipole  field 


,  „  ,  —Fo^fhu  sin  9  e^fc“r  7r/4f 

p(r,6,u!)~ - j== - ,  kwr 


V2 


oo 


nr 


(2.33) 


where  angle  6  is  defined  as  in  Figure  2.3.  Therefore,  to  better  visualize  the  emergence 
of  the  dipole  we  plot  the  directivity  (again,  oc  \p(r,  9,  cn)|2)  of  the  sound  normalized  with 
respect  to  its  value  at  0  —  |,  clipping  the  beams  in  the  fore  and  aft  directions,  as  in 
Figure  2.4. 


Figure  2.4.  Directivity  of  the  jet-excited  sound  for  the  general  case  of  an  air-filled 
spherical  cavity  at  STP  determined  by  Eq.  (2.32)  when  kwa  =  0.1,  1,  10,  100.  Each 
curve  is  normalized  with  respect  to  its  value  at  9  —  |. 


It  is  also  of  interest  to  compare  these  predictions  with  those  for  two  of  the  special 
cases  discussed  in  §2.1.3  and  §2.1.4.  In  §2.1.4  the  effect  on  sound  generation  of  wave 
incoherence  produced  by  a  randomly  irregular  cavity  interface  is  modeled  by  the  formal 
neglect  of  pressure  waves  incident  on  the  interface  from  within  the  cavity.  In  this  case 
\p(r,  d,cn)|2  in  the  far  held  is  given  by  equation  (2.32)  with  Z2n,2m  replaced  by  Z2npm 
defined  in  equation  (2.16).  The  corresponding  directivities  (in  Fig.  2.5)  are  seen  to  be 
hardly  changed  from  those  predicted  from  equation  (2.32)  for  the  full  solution,  indicating 
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that  the  directivity  of  high  frequency  resonances  are  essentially  unaffected  by  ‘wrinkling’ 
of  the  cavity  interface. 


Figure  2.5.  Directivity  of  the  jet-excited  sound  when  the  correction  of  Eq.  (2.16)  is 
introduced  into  (2.32)  to  account  for  incoherent  reflection  of  waves  from  the  interface,  for 
an  air-filled  spherical  cavity  at  STP  when  kwa  =  0.1,  1,  10,  100.  Each  curve  is  normalized 
with  respect  to  its  value  at  9  —  |. 


Second,  in  case  (ii)  of  §2.1.3  pressure  variations  within  the  cavity  are  absent  (a  ‘vac¬ 
uous’  cavity  interior).  It  now  follows  from  equation  (2.17)  that 


P(r,  0,  u) 

2  1 

\  ^  AnAmP 2n(cOS$).f) 2ttt,(cOS$) 

FQj  a 

~  (kwr)2 

h^(kwa)h^(kwa) 

(2.34) 


Again,  the  typical  field  shapes  plotted  in  Fig.2.6  are  effectively  identical  with  the  corre¬ 
sponding  plots  for  the  full  solution. 

This  agreement  occurs  because,  except  at  the  resonance  frequency,  Z2n/ (Z2n  —  1)  ~ 
Z2n/{Z2n-i)  ~  1  when  pw  p0 ■  This  conclusion  has  been  confirmed  (and  the  accuracy 
of  the  numerical  procedure  validated)  by  recalculating  the  directivities  for  the  fully  com- 
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1 


Figure  2.6.  Directivity  of  the  jet-excited  sound  in  case  (ii)  of  §2.1.3  (a  ‘vacuous’  cavity 
interior)  determined  by  Eq.  (2.34)  when  kwa  =  0.1,  1,  10,  100.  Each  curve  is  normalized 
with  respect  to  its  value  at  9  —  |. 

pressible  cavity  for  the  case  of  §2.1.4  of  incoherent  interface  scattering  in  cases  where 
Po  ~  pw  (but  with  unchanged  values  of  cQ  and  cw).  The  resulting  directivities  are  con¬ 
sistently  found  to  be  different  from  each  other,  and  also  different  from  the  results  in 
Figs. 2. 4  -  2.6  for  an  air  hlled  cavity. 

2.4  Conclusions 

The  low  frequency  sound  radiated  by  the  spherical  model  of  a  supercavity  is  dominated 
by  the  volumetric  monopole  modes  excited  by  the  jet.  The  amplitudes  of  motions  near 
the  cavity  are  controlled  by  damping  produced  by  the  sound  radiation,  by  viscous  and 
thermal  losses  in  the  cavity,  and  by  the  steady  hydrodynamic  flow  of  gas  from  the  ‘rear 
end’  of  the  cavity.  For  all  frequencies  of  practical  interest,  the  continuity  of  pressure 
across  the  cavity  interface  and  within  the  exhaust  flow  ensures  that  the  hydrodynamic 
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damping  is  negligible,  and  therefore  that  the  acoustics  of  the  system  are  effectively  un¬ 
affected  by  the  gas  exhaust  flow. 

The  radiation  from  the  cavity  is  omnidirectional  at  low  frequencies,  comparable  to  the 
Minnaert  frequency  of  the  dominant  monopole.  As  the  frequency  increases  the  radiation 
directivity  develops  a  principal  side-lobe  equivalent  to  that  produced  by  a  dipole  acoustic 
source  with  axis  along  the  radial  direction  of  the  gas  jet  impinging  on  the  interface.  Our 
numerical  results  for  the  spherical  cavity  with  a  jet  impinging  at  Q  —  \  to  the  nominal 
direction  of  the  mean  water  flow  indicate  that  this  dipole  radiation  (of  directivity  ~  sin2  6) 
is  accompanied  by  narrow  beams  radiated  in  the  fore  and  aft  directions  (6  =  0,  7r).  These 
beams  are  an  artifact  of  the  spherical  geometry  that  is  likely  to  be  absent  in  more  realistic 
cavity  geometries. 

The  theoretical  results  for  the  spherical  cavity  indicate  that  the  side-line  dipole  is 
well  formed  when  kwa  >  10.  For  the  more  general  cavity  this  suggests  that  the  dipole 
becomes  prominent  at  frequencies  exceeding  about  1.5cw/7?.Hz,  where  1Z  is  the  mean 
radius  of  curvature  of  the  interface  at  the  point  of  impact  of  the  jet.  At  such  frequencies, 
the  self-noise  generated  by  the  gas  jet  might  be  expected  to  have  a  negligible  influence 
on  control  systems  situated  forward  of  the  cavitator. 
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Chapter  3 


Analytical  model  for  jet  impingement 
on  a  plane  interface 


This  chapter  describes  the  sound  generated  by  a  circular  jet  of  small  cross-section  im¬ 
pinging  on  a  planar  interface.  We  consider  the  scale-model  vehicle  depicted  in  Figure  1.1. 
Ventilating  gas  is  injected  into  the  cavity  via  an  axisymmetric  ring  of  twenty  radiating 
jets  which  impinge  on  the  cavity  wall  at  normal  incidence.  Each  jet  impact  is  discreet 
and  will  be  considered  as  an  independent  source  of  sound.  Given  that  each  jet  is  identical 
and  the  injection  ring,  vehicle,  and  supercavity  are  all  axisymmetric,  it  is  therefore  only 
necessary  to  compute  the  sound  generated  by  the  impingement  of  a  single  jet. 

Next  we  consider  that  the  estimated  radius  of  curvature  of  the  supercavity  is  more 
than  three  times  the  diameter  of  the  jet  impact  area,  D.  Furthermore,  at  the  highest 
frequencies  u  of  interest,  kwD  ~  1,  where  acoustic  wavenumber  kw  =  co/cw  and  cw  is  the 
speed  of  sound  in  water,  so  we  will  approximate  the  source  as  being  acoustically  compact 
(see  Appendix  B).  We  can  therefore  approximate  the  site  of  the  jet  impact  on  the  cavity 
wall  as  locally  plane,  and  the  jet  to  be  of  infinitesimal  thickness.  As  with  the  sphere 
problem  of  the  previous  chapter,  these  simplifications  eliminate  difficulties  associated 
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with  a  more  realistic  cavity  geometry.  The  predicted  pressure  field  and  its  expression  in 
terms  of  the  frequency  spectrum  of  the  force  applied  by  the  jet  on  the  wall  is  determined. 


3.1  Pressure  fluctuations  in  the  water 


interface 


Figure  3.1.  Problem  of  a  jet  of  infinitesimal  cross-section  impinging  on  a  gas- water 
interface 

Let  the  mean  position  of  the  gas-water  interface  coincide  with  the  plane  X3  =  0  of 
the  rectangular  coordinate  system  x  =  (xi,X2,£3),  with  water  of  density  pw  and  sound 
speed  cw  occupying  the  region  X3  >  0  (Fig.  3.1).  The  Mach  number  of  the  mean  flow  of 
water  over  the  interface  is  assumed  to  be  sufficiently  small  that  the  convection  of  sound 
can  be  ignored.  Small  amplitude  fluctuations  in  the  pressure  p(x,  t )  (where  t  denotes 
time)  may  therefore  be  taken  to  satisfy 

(4^-v2)p  =  °-  >0-  P'1' 

provided  that  turbulence  sources  of  ‘noise’  within  the  flow  are  ignored  [18]. 

A  circular  gas  jet  which  is  coaxial  with  the  X3  axis  is  incident  normally  on  the  interface 
from  within  the  cavity  X3  <  0.  The  jet  exhausts  at  mean  speed  UQ  from  a  circular  nozzle 
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of  radius  Ra  whose  exit  plane  is  at  x3  =  —  £  (where  R0  £).  The  unsteady  force  F(t),  say, 
exerted  by  the  jet  on  the  interface  is  distributed  over  a  circularly  symmetric  region  of  the 
XiX2-plane  centered  on  the  origin  of  nominal  effective  radius  R  -C  t.  This  force  excites 
unsteady  motions  and  pressure  fluctuations  within  the  water.  When  attention  is  confined 
to  the  unsteady  motions  occurring  at  distances  |x|  R  and  when,  in  addition,  the 
relevant  acoustic  wavelengths  are  much  larger  than  R.  it  is  permissible  for  the  purposes 
of  calculation  to  assume  that  F(t)  is  concentrated  at  the  origin,  and  corresponds  to  a 
singular  distribution  of  applied  surface  pressure  equal  to  F(t)S(xi)S(x2)  ■  Furthermore, 
the  very  large  differences  in  the  mean  gas  and  water  densities  implies  that  to  an  excellent 
approximation  pressure  fluctuations  occurring  elsewhere  on  the  interface  can  be  ignored 
(for  further  discussion  of  this  approximation  see  Chapter  2).  The  solution  of  equation 
(3.1)  is  therefore  required  in  x3  >  0  subject  to  the  condition 

p  =  F(t)S(x1)d(x2)  on  x3  =  0.  (3.2) 


To  determine  p  we  introduce  a  Green’s  function  G(x,  y,t  —  r)  defined  by 


G(x,y,f  -  r) 


47t|x  —  y| 
1 


47t|x  —  y  | 


5\  t  —  T  — 


5  \  t  —  r  — 


(3.3) 


where  y  =  (yi,y2,  ~y3)  is  the  ‘image’  of  the  point  y  in  the  plane  of  the  interface.  The 
Green’s  function  vanishes  when  either  of  x3,  y3  lies  on  the  interface,  and  satisfies 


y  Cq  dt2 


G  =  6(x-  y)S(t  -  r), 


x3,  2/3  >  0, 


(3.4) 


where  G  =  0  for  t  <  r  (see  Appendix  B  for  a  brief  discussion  of  the  derivation  and  use 
of  Green’s  functions). 
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The  application  of  Green’s  theorem  and  the  radiation  condition  (see  Appendices  A 
and  B)  to  Eqs.  (3.1),  (3.4)  now  permits  the  pressure  in  x3  >  0  to  be  expressed  in  the 
form  [2,  6,  10,  29] 

P(x,t)  =  j>  fp(y,T)^(x,y,£-T) 

-G(x,y,t-T)J^(y,T)^  dy1dy2dr  (3.5) 

where  the  integration  is  over  all  values  of  the  retarded  time  — oo  <  r  <  +oo,  and  the 
surface  integral  is  over  the  cavity  interface  S:  y 3  =  +0. 

Substituting  into  the  integrand  from  Eqs. (3.2),  (3.3),  we  find  that  both  the  near  and 
far  field  pressure  fluctuations  in  the  water  are  given  without  further  approximation  by 
the  dipole  formula 


P(x,t) 


o  r  1 

8x3  27t|x 


F 


x3  >  0. 


(3.6) 


3.2  The  pressure  spectrum 

It  will  be  assumed  that  the  unsteady  component  of  the  jet  impact  force  F(t)  (the  net 
force  minus  the  mean  force)  is  stationary  random  with  frequency  spectrum 

=  ^~  [°°  ( F(t)F(t  +  t))  eiuJT  dr  (3.7) 

2tt  J- 00 

where  the  angle  brackets  (  )  denote  an  ensemble  or  time  average.  Then  (3.6)  implies 
that  the  corresponding  pressure  spectrum  $pp(x,  u)  in  the  water  is  given  by 
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2 


<f cos2  9 


x3  >  0, 


(3.8) 


47 r2  Ixl4 


where  cos  6  =  0:3/ |x|  determines  the  dipolar  directional  characteristics  of  both  the  near 
and  far  pressure  fields  in  the  water,  and  kw  —  uj/cw  is  the  acoustic  wavenumber,  and 


/ 


(p(x,  t)p(x,  t  +  t))  ei0Jr  dr. 


(3.9) 


—  OO 


Equation  3.8  determines  the  pressure  spectrum  at  x  in  the  water  in  terms  of  the 
spectrum  of  the  impact  force  of  the  jet.  In  the  next  chapter,  this  formula  will  be  used 
to  approximate  the  spectrum  of  the  water-borne  sound  generated  by  the  experimental 
vehicle  (Figure  1.1)  by  determination  of  the  impact  force  spectrum  exerted  by  one 
of  the  vehicle’s  jets. 
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Chapter  4 


Jet  impingement  experiment 


This  chapter  describes  a  canonical  experiment  performed  to  determine  the  frequency 
spectrum  &ff  of  the  interaction  force  for  a  single  jet  impinging  on  an  interface.  We 
again  consider  the  vehicle  depicted  in  Figure  1.1.  The  ventilating  gas  enters  the  cavity 
axisymmetrically  through  twenty  radially  orientated  nozzles  equally  spaced  on  a  gas 
injector  ring  close  to  the  cavitator.  Each  jet  impinges  on  the  cavity  interface  at  close  to 
normal  incidence.  In  the  impact  region,  turbulence  and  other  water-borne  disturbances 
such  as  bubbles  are  rapidly  convected  away  over  a  pressure  release  interface  by  the 
mean  flow  of  water  outside  the  experimental  cavity  and  exert  a  negligible  surface  force 
and,  therefore,  are  expected  to  have  no  effective  influence  on  acoustic  generation  and 
transmission.  The  wakes  behind  the  ridges  in  the  cavity  created  by  the  impinging  jets 
(see  Fig.  1.3)  are  not  acoustic  sources  because  they  convect  along  with  the  mean  flow 
velocity  and  exert  no  pressure  on  the  water.  Forces  exerted  by  turbulence  and  bubbles 
convecting  over  these  wakes  are  identically  canceled  by  equal  and  opposite  images  inside 
the  cavity  (see  Appendix  B).  Therefore,  the  dominant  acoustic  source  associated  with 
ventilation  jet  impingement  is  the  unsteady  fluctuations  of  the  cavity  wall  in  the  jet 
impact  region. 
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The  sound  produced  by  these  interactions  is  determined  principally  by  the  unsteady 
forces  exerted  on  the  water  at  the  interface  by  the  impinging  jets.  In  practice  the  ven¬ 
tilating  jets  are  fully  turbulent  on  impact  with  the  interface,  and  the  properties  of  the 
sound  depend  on  the  spectral  content  of  the  surface  force.  Typical  gas  speeds  on  im¬ 
pact  do  not  normally  exceed  ~  60  m/s.  In  these  circumstances  the  very  large  difference 
in  the  densities  of  the  gas  and  water  permit  the  interface  to  be  regarded  as  rigid  in  a 
first  approximation.  The  unsteady  impact  force  of  the  jet  may  then  be  identified  with 
that  obtained  when  a  jet  is  allowed  to  impinge  on  a  rigid  plate.  To  exploit  this  analogy 
one  of  the  ventilating  jets  of  the  test  vehicle  was  directed  at  normal  incidence  onto  the 
plate-like  measurement  surface  of  a  rigidly-fixed  force  transducer.  Though  made  slightly 
awkward  by  the  vehicle  body  and  mounting  system,  this  method  allowed  us  to  capture 
the  properties  and  frequency  dependence  of  the  vehicle  jet  due  to  its  unique  orifice  and 
injection  system,  details  which  may  be  lost  when  using  a  simpler,  stand-alone  jet  nozzle. 
The  spectrum  of  the  exerted  force  was  measured  for  a  range  of  gas  injection  flow  rates 
and  jet  lengths.  The  results  of  these  tests  were  used  with  the  formulae  determined  in 
Chapter  3  to  make  predictions  for  the  sound  generated  by  jet  impingement  on  a  planar 
interface. 

As  discussed  in  Section  1.2,  a  literature  survey  revealed  that  investigations  of  the 
acoustics  of  jet  impingement  [28,  34,  40]  have  focused  on  impingement  on  a  plate  and  on 
acoustic  sources  such  as  shed  vorticity  or  turbulence;  no  work  appears  to  have  been  done 
on  the  acoustics  of  a  jet  impinging  on  an  unsteady,  free  gas- water  interface.  Furthermore, 
no  examples  of  the  direct  measurement  of  the  unsteady  force  exerted  by  an  impinging 
jet  were  found.  Thus  the  few  existing  studies  of  relevance  are  those  which  provide 
information  about  the  exerted  pressure.  Strong  et  al.  [45]  describe  an  experiment  in 
which  a  jet  was  directed  at  a  flat  plate  at  various  angles  of  incidence.  Two  surface 
pinhole  microphones  were  used  to  measure  the  correlation  spectrum  of  the  fluctuating 
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surface  pressure  on  the  plate.  Their  measurements  for  normal  incidence  are  compared  in 
this  chapter  with  our  experimental  results  for  the  surface  force  spectrum  on  the  flat  force 
transducer.  We  also  compare  predictions  of  the  impingement-generated  sound  using  the 
data  from  [45]  with  our  predictions  for  the  supercavity. 

The  first  section  (4.1)  describes  the  supercavitating  vehicle  experimental  setup.  Sec¬ 
tion  4.2  recounts  the  experimental  method  and  executed  trials.  Data  analysis  and  un¬ 
certainty  in  the  results  are  discussed  in  Section  4.3.  Results  for  the  exerted  force  spectra 
are  presented  in  Section  4.4.  Using  these  results,  predictions  are  made  in  Section  4.5 
for  the  sound  generated  by  jet  impingement  on  a  planar  interface  and  compared  with 
corresponding  predictions  based  on  the  results  of  Strong  et  al.  [45]  in  Section  4.6.  Con¬ 
clusions  are  then  inferred  for  the  sound  generated  by  jet  impingement  on  a  more  realistic, 
cylindrical  cavity  interface. 

4.1  Experimental  setup 

Testing  was  conducted  at  the  Garfield  Thomas  Water  Tunnel  facility  at  ARL  using  a 
single  ventilating  jet  of  a  76  mm  body-diameter  experimental  supercavitating  vehicle  like 
that  depicted  in  Figures  1.1  and  1.2.  In  normal  operation  the  vehicle  is  rigidly  mounted 
in  an  ARL  0.305  m  water  tunnel  test  section  (see  Appendix  C).  Gas  is  injected  into  the 
cavity  axisymmetrically  via  a  ring  of  20  identical,  radially  directed  jets.  The  diameter 
of  the  exit  orifice  for  each  jet  is  1.6  mm,  and  the  exit  planes  are  evenly  distributed  in  a 
ring  of  radius  25  mm  centered  on  the  axis  of  symmetry  of  the  vehicle. 

Gas  enters  the  system  by  way  of  a  line  passing  through  the  vehicle  mounting  strut 
(Fig.  4.1).  A  line  from  each  of  six  bottles  of  compressed  gas  is  connected  to  a  manifold 
where  the  separate  supplies  are  combined  into  single  gas  line.  A  seventh  bottle  of  com¬ 
pressed  air  provides  a  constant  back  pressure  to  this  gas  line  at  the  manifold  so  that  a 
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constant  flow  rate  of  gas  to  the  vehicle  can  be  maintained  as  the  supply  bottles  empty.  A 
ball  valve  in  the  supply  line  allows  for  rapid  on/off  control  of  the  gas  flow  to  the  vehicle 
when  the  bottles  are  open  and  feeding  the  line.  The  gas  flow  rate  through  the  system 
is  monitored  by  a  flow  meter  upstream  of  the  ball  valve,  and  is  varied  by  adjustment  of 
the  ball  valve  aperture. 

back  pressure 
gas  bottle 


I 

I 

I 

manifold  ■ 


Figure  4.1.  Schematic  of  gas  injection  system 

The  model  vehicle,  without  the  cavitating  nose  section  which  was  not  required  for 
this  experiment,  was  secured  to  a  bench  top.  The  gas  injection  system  described  above 
was  connected  to  the  vehicle.  Four  of  the  jet  orifices  were  plugged  with  wax  held  in  place 
with  firmly  fixed  electrical  tape,  two  each  immediately  clockwise  and  anti- clockwise  of 
the  ‘test’  jet  (see  Fig  4.2).  Though  this  increased  the  net  flow  through  the  remaining 
open  orifices  by  approximately  25%  (relative  to  typical  running  conditions  with  all  jets 
operating),  it  was  important  to  do  this  in  order  to  isolate  a  single  jet  for  measurement. 
This  isolated  test  jet  was  directed  at  normal  incidence  onto  the  12.7  mm  diameter  mea¬ 
surement  surface  of  a  PCB  ICP  dynamic  force  transducer  (Fig.  4.3). 
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0  -  closed  orifice  0  -  test  jet  orifice 
0  -  open  orifice  []- transducer 


Figure  4.2.  Arrangement  of  open  and  closed  jet  orifices 


accelerometer 
test  jet  orifice 
yoke 

force  transducer 
plate 
vehicle 
translation  table 


Figure  4.3.  Experimental  setup.  The  vehicle’s  nose  section,  unnecessary  for  this  ex¬ 
periment,  has  been  removed.  In  this  photograph,  only  the  test  jet  is  exposed;  the  other 
fifteen  jets  were  subsequently  unsealed  as  in  Fig.  4.2  before  testing  was  conducted. 
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The  force  transducer  was  fixed  at  the  center  of  a  6.35nnn  thick  76.2  mm  square 
aluminum  plate  (Fig.  4.4).  The  plate  was  screwed  onto  a  127nnn  tall  heavy  aluminum 
yoke  made  massive  in  an  effort  to  reduce  structural  vibrations  induced  by  the  impacting 
jet  (Fig.  4.5).  The  yoke  was  fastened  on  a  1-axis,  0.25  mm  resolution  translation  table 
mounted  on  the  bench  top.  The  axis  of  the  translation  table  was  parallel  with  the  center 
axis  of  the  jet  allowing  for  adjustment  of  the  length  of  the  test  jet,  defined  to  be  the 
distance  between  the  jet  orifice  and  the  wetted  surface  of  the  transducer.  Superglue  was 
used  to  fix  an  accelerometer  to  the  plate,  in  close  proximity  with  the  transducer. 


Figure  4.4.  Plate  and  transducer  assembly.  The  force  transducer  was  threaded  into 
the  center  of  the  plate.  The  accelerometer  fixed  with  superglue. 


Figure  4.5.  Massive  yoke  used  for  mounting  of  plate  and  transducer  assembly 
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4.2  Methods 

The  PCB  dynamic  force  sensor  directly  measured  the  unsteady  force  exerted  by  the  test 
jet.  Pertinent  specifications  for  this  and  all  of  the  electronic  equipment  used  are  pre¬ 
sented  in  Appendix  D.  The  force  transducer  was  calibrated  using  a  set  of  known  masses. 
This  piezoelectric  transducer  cannot  measure  static  loads,  so  the  calibration  technique, 
similar  to  that  used  by  the  manufacturer,  consisted  of  securing  the  transducer  to  the 
benchtop  with  its  measurement  surface  horizontal  (i.e.  facing  upward),  resting  a  mass 
on  the  transducer,  and  measuring  the  output  voltage  when  the  mass  was  quickly  re¬ 
moved.  This  was  executed  three  times  each  for  masses  between  0.05  and  2  kg.  Assuming 
the  force  exerted  on  the  transducer  by  each  mass  was  its  weight,  mg,  these  data  and  the 
resulting  calibration  curve  are  shown  in  Figure  4.6.  The  calibration  curve  was  fitted  to 
the  data  using  the  method  of  linear-least-squares. 


Figure  4.6.  Calibration  data  for  unsteady  force  transducer 

The  accelerometer  was  used  to  measure  the  structural  vibrations  of  the  plate  gener¬ 
ated  from  both  the  test  environment  and  the  jet  impact.  Assuming  a  rigid  mounting  (the 


38 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


force  transducer  was  threaded  tightly  onto  the  plate),  this  vibration  is  taken  to  be  that 
experienced  by  the  force  transducer  housing.  The  accelerometer  was  calibrated  using  a 
PCB  394C06  hand  held  1-g  shaker.  The  shaker  is  an  oscillating  piston  which  accelerates 
the  transducer  at  1  g.  The  accelerometer  was  fixed  to  the  shaker  and  the  amplitude  of 
its  output  voltage  recorded.  This  was  repeated  three  times  and  the  results  averaged, 
yielding  a  sensitivity  of  14.566  mV/g. 

The  signals  from  the  force  transducer  and  accelerometer  were  fed  through  a  PCB 
480B21  signal  conditioner,  applying  a  gain  of  10  to  both  the  force  sensor  signal  and 
the  accelerometer  signal.  The  signals  were  then  passed  to  a  DSP  SigLab  50-21  data 
acquisition  system.  This  system  was  connected  to  a  laptop  computer  with  accompanying 
SigLab  software  used  to  record  and  process  each  signal.  The  software  applied  a  fast- 
Fourier  transform  (FFT)  of  8192  points  to  64  ms  samples  of  the  time  series  of  each  signal 
to  produce  an  autocorrelation  spectrum  with  units  of  root-mean-square  (RMS)  voltage 
squared  in  15.6  Hz  bins  up  to  50  kHz  (see  Bendat  and  Piersol  [3]  for  a  discussion  of  the 
theory  and  application  of  spectral  functions).  An  ensemble  of  60  of  these  spectra  were 
averaged  to  produce  a  final  resultant  spectrum.  Cross-correlation  spectra  between  the 
force  sensor  and  accelerometer  signals  were  also  generated  using  the  same  transform  and 
averaging  parameters. 

Measurements  were  made  for  a  variety  of  jet  lengths  i  and  gas  injection  volume  flow 
rates  Q.  A  complete  list  of  trial  parameters  is  given  in  Table  4.1.  Typical  experimental 
injection  rates  for  the  model  test  vehicle  are  0.0047,  0.0071,  and  0.0094  m3/s,  producing 
jet  lengths  of  approximately  35nnn.  A  range  of  lengths  around  this  median  value  were 
tested  to  study  the  effect  of  jet  length  on  the  spectrum  of  the  applied  force  and  the 
produced  sound.  For  the  shorter  of  these  jet  lengths,  15  and  20  mm,  the  entire  impact 
area  of  the  jet  lies  within  the  face  of  the  force  transducer,  assuming  a  jet  spread  half-angle 
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of  12.5°  [29]  and  finding  the  diameter  D  of  the  jet  cross-section 

D  —  2(i?0  +  f?tan(12.5°)),  (4.1) 

in  terms  of  the  radius  of  the  jet  orifice  Ra  and  jet  length  i  (see  Figure  4.7).  Longer 
jet  lengths  have  cross-sections  on  impact  larger  than  the  force  transducer  measurement 
surface.  However,  the  bulk  of  the  force  exerted  by  a  jet  is  produced  by  the  flow  in  the 
vicinity  of  the  jet  centerline,  inferred  by  the  variation  of  longitudinal  velocity  with  radial 
distance  from  the  jet’s  central  axis.  If  the  jet’s  velocity  drops  off  rapidly  above  0.4a  from 
the  centerline  [29],  we  estimate  from  (4.1)  that  the  force  transducer  is  large  enough  to 
capture  the  bulk  of  the  force  exerted  by  a  jet  up  to  64  mm  in  length.  It  was  therefore 
assumed  that  for  our  jets  of  length  >  20  mm,  the  force  measured  by  the  transducer  is  a 
sufficient  approximation  to  the  force  exerted  over  the  entire  impact  area.  Comparison  of 
the  measured  force  spectra  for  the  longer  jet  lengths  with  those  for  jet  lengths  of  15  and 
20  mm  will  demonstrate  the  appropriateness  of  this  assertion  and  the  magnitude  of  any 
introduced  error. 


i  (mm) 

Q  (m3/s)  xlO4 

15 

47,  71,  94 

20 

47,  71,  94 

25 

47,  71,  94 

30 

47,  71,  94 

35 

47,  71,  94 

40 

47,  71,  94 

45 

47,  71,  94 

Table  4.1.  Schedule  of  experimental  trials 

Each  trial  was  repeated  three  times.  In  addition,  the  centerline  jet  exit  velocity 
UQ  at  the  jet  orifice  was  measured  for  a  jet  injection  rate  of  Q  =  0.0047m/s  using  a 
TSI  VelociCalc  probe  air  velocity  meter.  Direct  measurement  of  the  exit  velocity  for 
the  higher  injection  rates  was  not  possible  because  of  the  limitations  imposed  by  the 
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Figure  4.7.  Estimated  shape  of  a  fully  turbulent  jet 

velocimeter’s  usable  range;  these  values  will  be  interpolated  from  the  lower  injection  rate 
measurement  using  fluid  dynamic  theory.  To  make  the  measurement,  the  probe  tip  was 
inserted  into  the  center  of  the  jet  approximately  1  mm  downstream  of  the  jet  orifice. 
The  VelociCalc  is  calibrated  annually  by  the  manufacturer,  and  the  measurement  was 
repeated  three  times  to  account  for  statistical  variability.  The  device  is  also  equipped 
with  a  temperature  probe  which  was  used  to  measure  the  ambient  air  temperature  and 
that  of  the  jet  at  the  jet  orifice  employing  the  same  technique. 

4.3  Data  analysis 

4.3.1  Post-processing 

The  determined  calibration  factor  for  each  transducer  was  applied  during  post-processing 
to  convert  the  measured  autocorrelation  spectra  from  units  of  RMS  voltage  squared  into 
units  of  squared  force  (N2)  or  acceleration  (g2)  and  the  cross-correlation  spectra  from 
units  of  RMS  voltage  squared  into  units  of  force x  acceleration,  N-g’s.  Applying  a  method 
similar  to  that  used  by  Young,  et  al.  [48],  the  spectra  were  then  combined  to  eliminate 
vibrational  contamination  in  the  measured  unsteady  force  spectra  (the  source  of  this 
contamination  is  discussed  later).  A  cross  correlation  spectrum  Xap{uS)  (where  angular 
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frequency  oj  =  2irf)  between  two  signals  a(t)  and  (3{t)  is  defined  by  [3] 

Xa0(uj)  =  [  {a(t)P(t  +  r))eUi)Tdr.  (4.2) 

The  coherence  between  the  signals  is  then  given  (in  dimensions  of  a2)  by  [3] 

Xap(u)  ■  Xpa(u)/App(u),  (4.3) 

where  App(uj)  =  Xpp(uj)  is  the  autocorrelation  spectrum  of  signal  /?.  Applying  this  to 
the  experiment,  vibrational  contamination  was  approximately  eliminated  by  subtraction 
from  the  measured  unsteady  force  spectrum  of  any  coherence  with  the  vibration  spectrum 
(measured  by  the  accelerometer),  according  to 

A ff,c(w)  =  Aff(lu)  —  XaF(uj)  ■  XFa{u)/ Aaa(u),  (4.4) 

where  AFF)C(uj)  is  the  ‘corrected’  unsteady  force  spectrum,  Aff(uj)  and  Aaa(co)  are  re¬ 
spectively  the  measured  unsteady  force  and  vibration  spectra,  and  XaF(u j)  and  XFa(ui) 
are  cross-correlation  spectra  between  the  accelerometer  and  force  transducer  signals. 

4.3.2  Uncertainty 

The  major  source  of  error  in  these  measurements  is  vibrational  contamination  of  the  force 
sensor  readings.  Two  sources  of  this  contamination  exist:  background  vibration  from  the 
experimental  environment  and  structural  vibrations  in  the  transducer  mounting  assembly 
induced  by  the  test  jet.  Transducer  readings  were  collected  with  the  jet  off  (i.e.  Q  =  0) 
to  record  any  environmental  noise,  as  well  as  electrical  noise  inherent  in  the  transducers 
or  electronics.  These  readings  were  consistently  10  to  40  dB  below  measurements  taken 
with  the  jet  on  except  at  the  lowest  (<  100  Hz)  and  highest  frequencies  (>  30kHz). 
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Examples  of  the  unsteady  force  and  vibration  spectra  recorded  for  a  jet  with  Q  =  0  and 
Q  =  0.0047  m3/s,  i  =  35  mm  are  given  in  Figures  4.8  and  4.9. 


10  100  1000  10000  100000 


frequency  (Hz) 


Figure  4.8.  Comparison  of  the  vibration  spectra  for  jets  of  flow  rates  Q  —  0  and 
Q  =  0.0047  m3/s  and  jet  length  i  =  35  mm.  The  spectra  are  given  in  units  of  dB  relative 
to  1  g2. 


Jet  induced  vibrations  were  partially  eliminated  by  application  of  the  coherence  cor¬ 
rection  formula  (4.4).  Figure  4.10  presents  an  example  of  the  unsteady  force  spectrum 
recorded  for  a  trial  at  Q  —  0.0047  m3/s  and  £  =  35  mm  before  and  after  application  of 
the  coherence  formula.  Noise  up  to  20dB  was  removed  at  frequencies  across  the  entire 
measured  bandwidth.  However,  it  appears  that  there  is  still  significant  contamination  of 
this  unsteady  force  spectrum  above  3000  Hz.  Comparing  corrected  spectra  for  varying 
jet  lengths  and  injection  rates  (see  an  example  in  Figure  4.11),  we  find  that  this  con¬ 
tamination  is  consistently  present  in  magnitude  and  frequency  dependence  in  all  of  the 
readings,  implying  that  the  contamination  is  due  to  structural  vibrations,  rather  than 
a  true  measurement  of  the  unsteady  force  exerted  by  the  jet.  We  will  therefore  choose 
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frequency  (Hz) 

Figure  4.9.  Comparison  of  the  unsteady  force  spectra  for  jets  of  flow  rates  Q  =  0  and 
Q  =  0.0047  m3/s  and  jet  length  t  =  35  mm.  The  spectra  are  given  in  units  of  dB  relative 
to  1  N2. 


to  estimate  the  unsteady  force  over  these  higher  frequencies  by  projecting  a  curve  which 
continues  the  slope  of  the  spectrum  at  lower  frequencies  (see  Figure  4.12). 

The  normalized  random  statistical  error  r\r  for  each  spectrum  can  be  estimated  as  [3] 

Vr  —  —j=  (4.5) 

Jn 


where  n  is  the  number  of  samples  averaged  in  the  ensemble.  As  given  in  Section  4.2, 
60  samples  were  averaged  in  each  ensemble  for  an  error  of  ±0.13.  Random  error  in  the 
readings  as  assessed  by  variation  between  duplicated  runs  (repeatability)  is  estimated  at 
±1.5 dB  (see  Fig.  4.13  for  an  example).  Errors  in  frequency  are  systematic,  dominated 
by  the  chosen  bin  width  of  the  Fourier  transform,  and  are  on  the  order  of  ±7.8  Hz.  The 
combined  statistical  and  bias  errors  in  the  jet  exit  velocity  measurements  produced  an 
error  of  ±1  m/s. 
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frequency  (Hz) 


Figure  4.10.  Comparison  of  the  corrected  and  uncorrected  unsteady  force  spectra 
(. Aff,c  and  App  of  Eq.  4.4,  respectively)  for  Q  =  0.0047  m3/s  and  i  =  35  mm.  The 
spectra  are  given  in  units  of  dB  relative  to  1  N2. 


Finally  we  consider  the  validity  of  the  force  spectra  measured  for  jet  lengths  >20  mm, 
for  which  the  cross-sectional  area  of  the  jet  at  impact  is  larger  than  the  force  transducer 
measurement  surface.  Figure  4.14  compares  the  corrected  force  spectra  for  trials  at 
each  jet  length  i  =  15,  20,  30,  and  40  mm  and  injection  rate  Q  =  0.0047m3/s.  The 
greatest  disparity  between  the  spectra  occurs  near  1000  Hz,  which  may  imply  that  these 
frequencies  are  dominant  at  the  outer  edge  of  the  jet’s  cross-section  and  they  are  not 
being  captured  by  the  transducer  at  longer  jet  lengths.  However,  across  the  rest  of  the 
measured  bandwidth  we  End  an  average  difference  of  approximately  ±2  dB  between  the 
spectra.  This  difference  is  of  the  same  order  as  the  statistical  error  between  repeated 
runs  (see  above);  we  therefore  conclude  that  the  force  spectrum  measurements  for  jets 
with  lengths  longer  than  20  mm  sufficiently  approximate  the  entire  force  exerted  by  the 
jet  across  these  frequencies.  Furthermore,  the  disparity  between  the  spectra  occurs  well 
below  the  frequencies  of  interest  for  the  self  noise  problem,  so  we  consider  it  acceptable 
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Figure  4.11.  Comparison  of  the  corrected  unsteady  force  spectra  for  jets  with  Q  = 
0.0047  and  0.0094  m3/s  and  i  =  15  and  35  mm.  The  spectra  are  given  in  units  of  dB 
relative  to  1  N2. 


to  use  these  measurements  to  make  predictions  for  the  sound  generated  by  the  longer 
jets. 
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frequency  (Hz) 


Figure  4.12.  Demonstrating  the  estimate  of  the  unsteady  force  spectrum  for  frequencies 
>3000  Hz.  The  given  spectrum  was  measured  for  a  jet  of  length  t  =  15  mm  and  injection 
rate  Q  =  0.0094  m3/s,  and  is  presented  in  units  of  dB  relative  to  1  N2. 


10  100  1000  10000  100000 


frequency  (Hz) 


Figure  4.13.  Comparison  of  corrected  force  spectra  for  multiple  trials  at  Q  =  0.0094 
m3/s  and  i  =  35  mm.  The  spectra  are  given  in  units  of  dB  relative  to  1  N2. 
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Figure  4.14.  Comparison  of  corrected  force  spectra  for  trials  at  Q  =  0.0047  m3/s  and 
l  =  15,  20,  30,  and  40  mm.  The  spectra  are  given  in  units  of  dB  relative  to  1  N2. 
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4.4  Results 


Figure  4.15.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s 
and  jet  lengths  £  =  15,  25,  35,  and  45  mm,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


The  corrected  unsteady  force  spectra  are  taken  to  be  the  spectra  of  the  impact  force 
for  a  single  jet  on  the  supercavity’s  gas-water  interface  (i.e.  =  Aff,c{ u;)).  Figures 

4.15-4.17  depict  the  measured  normalized  impact  force  spectrum 


10  x  log10 


(  ®ff(w)U0\ 

V  F20D0  ) 


as  a  function  of  non-dimensional  frequency  fD0/U0  for  the  gas  flow  rates  of  0.0047, 
0.0071,  and  0.0094  m3/s,  respectively,  where  Fa  is  the  mean  force  exerted  by  the  jet,  Da 
is  the  diameter  of  the  jet  orifice,  and  UQ  is  the  mean  centerline  jet  velocity  at  the  orifice. 
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Figure  4.16.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s 
and  jet  lengths  £  =  15,  25,  35,  and  45  mm,  normalized  to  mean  force  Fa  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


F0  is  defined  by  the  momentum  balance  formula 


F0 


(4.6) 


where  U  is  the  mean  jet  velocity  at  impact,  D  is  the  diameter  of  the  jet  impact  area, 
and  pa  is  the  mean  density  of  the  jet.  The  density  pa  is  approximated  as  that  for  air  at 
STP,  1.2  kg/m3,  based  on  temperature  readings  of  the  jet  which  measured  only  about  3° 
lower  than  ambient.  D  is  calculated  for  each  jet  length  £  from  D0  =  2 R0  using  Eq.  (4.1). 
As  described  in  §4.2,  the  jet  exit  centerline  velocity  UQ  was  measured  for  the  jets  of  Q 
=  0.0047  m/s  with  a  probe  velocimeter.  The  value  of  U0  was  determined  for  higher  flow 
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Figure  4.17.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s 
and  jet  lengths  £  =  15,  25,  35,  and  45  mm,  normalized  to  mean  force  Fa  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


rates  by  assuming  a  linear  relationship  between  flow  rate  and  velocity,  i.e. 

U2  =  (4.7) 

VI 

where  U\  and  U2  are  the  exit  velocities  which  correspond  to  injection  flow  rates  Q\  and 
Q2,  respectively.  This  relationship  is  derived  using  the  definition  of  the  volume  flow  rate 
out  of  the  jet  orifice  Q0  =  U0(tt / F)D20. 

The  volume  flow  rate  Q'  =  U(ir/4)D2  through  any  cross-section  of  the  jet  at  length 
£  is  given  by  [29] 

Q'  =  PQo-^-  (4.8) 

1\0 

For  a  circular  jet  orifice,  the  empirical  value  of  j3  =  1.5  [29].  Applying  the  definitions  of 
Q' ,  Q0  to  (4.8),  the  mean  jet  impact  velocity  for  a  jet  of  length  £  is  interpolated  from  the 
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exit  velocity  using  [29] 


(4.9) 


This  self-similarity  relation  breaks  down  near  the  jet  orifice.  For  the  shortest  jet  length  t 
=  15  mm,  the  impact  velocity  is  unlikely  to  differ  signihcanly  from  the  exit  velocity  [16]; 
we  will  therefore  take  U  ~  UQ.  Values  of  D  and  U  ( D0  and  UQ  at  the  jet  orifice,  £  =  0) 
for  each  case  are  tabulated  below. 


i  (mm) 

0 

15 

20 

25 

30 

35 

40 

45 

D0,  D  (mm) 

1.6 

8.3 

10.5 

12.7 

14.9 

17.1 

19.3 

21.6 

u0, 

u 

(m/s);  Q 

=  0.0047  m3/s 

27.1 

27.1 

23.7 

20.2 

17.6 

15.5 

13.9 

12.6 

Uo , 

u 

(m/s);  Q 

=  0.0071  m3/s 

40.7 

40.7 

35.6 

30.3 

26.3 

23.3 

20.9 

18.9 

Uo, 

u 

(m/s);  Q 

=  0.0094  m3/s 

54.2 

54.2 

47.5 

40.4 

35.2 

31.1 

27.8 

25.2 

Table  4.2.  Estimated  values  of  jet  diameter  D  and  exit  (£  =  0)  and  mean  impact 
velocities  UQ  and  U  for  each  trial  at  jet  length  i  and  injection  flow  rate  Q. 


In  all  cases  the  spectra  in  Figures  4.15-4.17  peak  near  fD0/U0  ~  0.01  then  drop  off  rapidly 
like  ( fD0/U0)~ 2'5.  For  the  purpose  of  presentation,  the  spectra  for  the  intermediate  jet 
lengths  of  20,  30,  and  40  mm  have  been  omitted;  the  normalized  unsteady  force  spectra 
of  all  of  the  collected  data  are  presented  in  Appendix  E. 

Figure  4.18  presents  the  unnormalized  measured  force  spectra  plotted  against  fre¬ 
quency  for  jets  with  test  parameters  Q  =  0.0047  and  0.0094  m3/s  and  i  =  15  and  35  mm. 
This  plot  elucidates  the  relationship  between  injection  rate  and  the  magnitude  of  the 
measured  force  spectrum  as  well  as  jet  length  and  magnitude.  Let  us  consider  that  the 
unsteady  component  F  of  the  force  exerted  by  the  jet  scales  with  the  mean  component 
of  the  exerted  force  F0  and  the  unsteady  component  of  the  jet’s  impact  velocity  u  scales 
with  the  mean  component  of  the  impact  velocity  U.  As  &ff  oc  F2  oc  u4,  then  a  dou¬ 
bling  of  the  injection  flow  rate  (and  e.g.  the  mean  impact  velocity  U)  results  in  four-fold 
increase  in  the  mean  force,  for  a  sixteen- fold  or  12dB  increase  in  F 2  oc  &ff-  Our  data 
approximate  this  well,  with  an  estimated  10-15  dB  increase  in  the  unsteady  force  spec- 
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tram  with  the  doubling  of  the  injection  flow  rate.  There  is  also  only  a  2-3  dB  difference 
between  the  unsteady  force  spectra  for  jets  of  differing  length;  this  is  attributable  to 
statistical  variability  between  runs. 


frequency  (Hz) 


Figure  4.18.  Measured  force  spectra  for  jets  with  test  parameters  Q  =  0.0047  and 
0.0094  m3/s  and  £  =  15  and  35  mm.  The  spectra  are  given  in  units  of  dB  relative  to  1 


N2. 


4.5  Predicted  acoustic  power  spectra 


These  impact  force  spectra  can  be  used  with  the  results  of  Chapter  3  to  predict  the 
spectrum  of  the  sound  power  generated  by  jet  impingement  on  a  planar  interface.  From 
the  definition  of  sound  power  (see  Appendix  A),  the  radiated  sound  power  spectrum 
fl(a;),  the  power  spectrum  for  the  sound  which  transmits  to  the  far  held,  is  given  by 


nH=  r 

J  0  JO  PwC-w 


(4.10) 
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where  the  integration  is  over  the  surface  of  a  large  hemisphere  in  the  water  of  radius  |x 


centered  on  the  centroid  of  the  gas  interface  impact  region.  Note  that  in  (3.8)  the  terms 
in  parenthesis  represent  respectively  the  near  and  far  pressure  fields,  and  the  near  field 
term  vanishes  as  |x|  — >  oo,  yielding 


(4.11) 


A  convenient  nondimensional  form  is 


pwIl(uj)  _  ^FF(u)cwn(kwD)2 


(4.12) 


M2p\U2D3  F2D  96 


where  M  is  the  jet  Mach  number  U/cw  relative  to  the  speed  of  sound  in  water.  The 


variation  of 


with  reduced  frequency  kwD  is  plotted  in  Figures  4.19-4.22  for  the  measured  cases  i  = 
15,  25,  35,  and  45  mm  at  0.0047,  0.0071,  and  0.0094  m3/s.  Predictions  of  the  radiated 
sound  power  spectrum  based  on  all  of  the  measured  force  spectra  are  given  in  Appendix 


E. 


The  predictions  have  collapsed  to  a  common  form  across  all  frequencies,  peaking  in 


0.005  <  kwD  <  0.05  and  dropping  off  like  kwD  0  5.  This  implies  that  the  dominant 


sound  is  in  a  frequency  range  where  the  wavelength  A  =  2n/kw  is  much  larger  than  the 
characteristic  length  D,  validating  our  hypothesis  of  acoustic  compactness  required  by 
the  plane  interface  approximation.  The  relative  differences  between  runs  of  varying  flow 
lengths  are  attributable  to  statistical  variation  for  each  run  and  the  differing  magnitudes 
of  the  normalizing  quantities.  As  II  oc  &FF,  Figure  4.18  and  its  discussion  at  the  end  of 
§4.4  indicate  that  there  is  approximately  a  10-15  dB  increase  in  radiated  sound  power 
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Figure  4.19.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  £  =  15  mm  versus  reduced  frequency  kw/D. 

with  a  doubling  of  gas  injection  rate. 

4.6  Comparison  with  Strong  et  al. 

In  order  to  have  some  independent  support  for  these  results  we  look  to  the  work  of  Strong 
et  al.  [45].  In  their  experiment  an  air  jet  with  nozzle  diameter  102  mm  was  directed  at  a 
heavy  plate.  Two  pinhole  microphones  flush  mounted  in  the  plate,  one  at  the  centroid  of 
the  impact  area  and  the  second  at  various  radial  distances  from  this  centroid,  were  used 
to  measure  the  correlation  of  the  unsteady  surface  pressure  on  the  plate  for  frequencies 
between  40  and  10,000  Hz.  Their  results  include  the  frequency  spectral  density  of  these 
pressure  fluctuations  $p(u;)  (in  their  nomenclature)  on  the  plate  for  a  jet  at  normal 
incidence.  This  is  a  measurement  of  the  unsteady  pressure  exerted  on  the  plate  at  the 
centroid  of  the  impact  area,  not  a  measurement  of  the  pressure  spectrum  of  the  sound 
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Figure  4.20.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  d  =  25  mm  versus  reduced  frequency  kw/D. 

radiating  from  the  jet  impact  region,  which  we  have  previously  called  <f>Pp(c<;).  These 
data  are  reproduced  in  Figure  4.23,  normalized  by  mean  jet  density  p^.,  jet  exit  velocity 
U0,  and  jet  orifice  diameter  Da  and  plotted  against  the  dimensionless  frequency  u)D0/U0. 

The  pressure  spectrum  can  be  used  to  make  a  rough  estimate  of  the  spectrum  of 
the  force  exerted  on  the  plate  by  the  jet.  As  a  first  approximation  we  assume  that  all 
frequency  components  of  the  pressure  on  the  plate  are  fully  correlated  over  the  impact 
area  |H2,  i.e.  we  assume 

=  (4.13) 

This  cannot  be  correct,  because  the  correlation  area  of  the  surface  pressure  must  be 
smaller  than  jD2.  Therefore  the  approximation  (4.13)  will  provide  an  upper  bound 
for  the  acoustic  radiation  spectrum.  Again  using  Eqs.  (4.1)  and  (4.6)  to  define  the 
experiment’s  characteristic  dimensions  of  mean  impact  velocity  and  mean  force,  the 
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Figure  4.21.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  l  =  35  mm  versus  reduced  frequency  kw/D. 


normalized  force  spectrum  defined  by  (4.13)  becomes 


10  x  log10 


f  $ff(u)U0\ 

V  f*d0  ) 


=  10  X  log10 


V  p\uw0 ) 


which  is  plotted  in  Figure  4.24  versus  fD0/U0.  Also  plotted  are  measured  unsteady 
force  spectra  for  the  supercavity  jets;  because  the  experiment  in  Strong  et  al.  was 
conducted  with  a  jet  length  of  14R„,  we  have  chosen  the  most  similar  cases  of  l  =  15  mm 
(=  18.75.Ro)  and  Q  =  0.0047,  0.0071,  and  0.0094  m3/s  for  comparison.  As  with  the  ARL 
jets,  the  Strong  et  al.  estimated  unsteady  force  spectrum  drops  off  like  (fD0/U0)~2'5. 
The  apparent  shift  in  the  curve  is  likely  a  consequence  of  the  differing  normalization 
values  for  each  experiment,  while  the  similar  frequency  dependence  confirms  the  validity 
of  our  measurement  technique  and  estimation  for  higher  frequencies. 

The  estimated  force  spectrum  is  used  to  predict  the  sound  power  spectrum  generated 
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Figure  4.22.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  £  =  45  mm  versus  reduced  frequency  kw/D. 


when  the  jet  impinges  on  a  gas-water  interface,  given  by  Eq.  4.11.  Figure  4.25  presents 
the  normalized  sound  power  spectra 


10  x  log10 


(  Po nM  ^  JD 
\M2p\U2D 3) 


(plotted  against  dimensionless  frequency  kwD)  for  the  jet  of  Strong  et  al.  compared  with 
corresponding  predictions  based  on  our  measurements  using  a  supercavitating  vehicle 
jet.  The  Strong  et  al.  spectrum  peaks  near  kwD  =  0.2  and  drops  off  like  (kwD)~0'5. 
The  similarity  of  the  frequency  dependence  between  the  predictions  for  the  Strong  et  al. 
and  ARL  jet  again  confirms  the  validity  of  our  chosen  measurement  technique  and  the 
generality  of  our  mathematical  solution  for  the  acoustic  pressure. 
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coD  0  /U  „ 


Figure  4.23.  Pressure  spectral  density  at  the  centroid  of  the  impact  area  for  the  jet 
measured  by  Strong  et  al.  [45]. 

4.7  Conclusions 

Observations  about  the  noise  predictions  for  a  plane  interface  can  be  used  to  postulate 
similar  characteristics  for  the  sound  generated  by  a  ventilation  jet  impinging  on  a  more 
realistic  cylindrical  interface  of  a  supercavity.  The  observed  small  differences  in  force 
spectra  associated  with  differences  in  the  jet  length  indicate  that  for  an  experimental 
supercavity,  where  an  accurate  determination  of  the  jet  length  is  not  easily  made,  a 
rough  estimate  of  £  to  within,  say,  5  mm,  is  sufficient  to  derive  an  adequate  prediction 
of  the  radiated  sound. 

The  measured  effect  of  increased  flow  rate  on  the  magnitude  of  the  impact  force 
spectrum  for  a  single  jet  on  a  plane  interface  will  likely  be  representative  of  the  multiple 
jet  case  for  the  supercavity.  For  the  planar  interface  Fig.  4.18  and  Eq.  4.11  indicate  that 
the  sound  power  increases  by  10-15  dB  when  the  injection  flow  rate  is  doubled.  This 
may  not  be  the  case  for  multiple  jets  because  of  possible  correlations  between  the  surface 
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fD0/U„ 


Figure  4.24.  Comparison  of  unsteady  force  spectra  for  ARL  jets  of  Q  =  0.0047,  0.0071, 
and  0.0094  m3/s  and  i  =  15  mm  and  the  jet  of  Strong  et  al.  [45]. 

force  fluctuations  produced  by  neighboring  jets.  In  any  event,  because  an  increase  in 
the  gas  injection  rate  causes  a  rapid  increase  in  sound  power  it  would  evidently  be  of 
considerable  interest  to  determine  the  minimum  gas  injection  rate  that  can  safely  sustain 
the  supercavity  of  an  underwater  vehicle  in  a  stable  form. 

ft  can  be  anticipated  that  the  directivity  of  sound  radiated  from  a  real  supercavity 
must  differ  in  detail  from  that  from  a  planar  interface  because  of  possible  interference 
between  neighboring  jets  and  the  fundamental  geometric  differences.  On  the  planar  in¬ 
terface  the  jet  produces  a  single  dipole  source  projecting  into  the  water,  normally  to 
the  interface,  while  the  supercavity  creates  an  axisymmetric  array  of  dipoles  projecting 
normally  off  of  the  cavity  into  the  water.  However,  the  far  field  magnitudes  and  fre¬ 
quency  dependences  of  the  generated  sound  should  be  similar.  The  sound  power  spectra 
predicted  for  the  plane  interface  indicate  that  most  of  the  radiated  sound  is  contained 
in  the  frequency  range  0.005  <  kwD  <  0.05.  This  confirms  the  assumption  that  the  jet 
impact  region  is  acoustically  compact,  and  therefore  that  the  local  geometry  of  the  im- 
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kwD 


Figure  4.25.  Comparison  of  predicted  sound  power  spectra  for  ARL  jets  of  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  £  =  15  mm  and  the  jet  of  Strong  et  al.  [45]. 

pact  region  is  probably  irrelevant.  The  power  spectrum  for  the  sound  radiated  from  the 
supercavity  would  be  expected  to  exhibit  a  frequency  dependence  that  is  very  similar  to 
that  predicted  for  the  plane  interface,  again  with  the  peak  acoustic  frequencies  between 
0.005<  kwD  <  0.05. 
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Chapter  5 


Self-noise  at  the  vehicle  nose 


In  this  chapter  we  predict  the  high-frequency  impingement  generated  noise  at  the  vehicle 
nose.  We  choose  to  model  the  experimental  ventilated  cavity  as  a  gaseous,  circular 
cylinder  seated  behind  a  solid  ellipsoidal  cavitating  nose  section  (see  Figure  5.1).  The 
ventilating  jets  exert  an  unsteady  pressure  on  the  cavity  interface  downstream  of  the 
cavitator.  This  constitutes  an  acoustic  source  of  dipole  type  which  generates  sound  in 
the  water.  By  finding  a  suitable  frequency-dependent  transfer  function  we  can  relate  the 
jet-exerted  pressure  to  the  corresponding  acoustic  pressure,  or  ‘self- noise,’  at  the  vehicle 
nose.  The  transfer  function,  computed  for  the  simplified  ellipsoidal  nose  and  cylindrical 
cavity  geometry,  takes  account  of  diffraction  of  high  frequency  waves  and  creeping  waves 
on  the  nose  and  the  surface  transition  from  the  solid  nose  section  to  the  gas-water  cavity 
interface. 

The  solution  for  the  acoustic  pressure  at  the  nose  in  terms  of  the  transfer  function 
is  given  in  Section  5.1.  The  transfer  function  is  determined  analytically  in  Section  5.2 
and  approximated  numerically  for  high  and  intermediate  frequencies  in  Section  5.3.  In 
Section  5.4,  the  spectrum  of  the  sound  at  the  nose  is  derived  in  terms  of  the  spectrum  of 
the  unsteady  force  exerted  by  a  ventilating  jet.  In  Section  5.5,  the  results  of  §5.3,  §5.4 
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are  combined  with  the  experimental  results  of  Chapter  4  to  predict  the  self-noise  at  the 
nose  generated  by  the  experimental  supercavitating  vehicle. 


5.1  Pressure  fluctuations  at  the  nose 


sound 


relative  mean  flow 
- >■ 


Figure  5.1.  Model  of  a  cylindrical  supercavity  with  an  ellipsoidal  cavitator. 


We  consider  the  undisturbed  gas-water  interface  (Sc  in  Figure  5.1)  to  be  a  circular 
cylinder  of  radius  r  =  a  coaxial  with  the  positive  x  axis,  taking  the  origin  at  the  midpoint 
O  of  the  rear,  plane  face  of  the  cavitator.  The  surface  Sn  of  the  cavitator  in  contact  with 
the  water  is  defined  in  terms  of  the  cylindrical  coordinates  of  Figure  5.1  by  the  ellipsoid 


r 


—b<x<  0,  b  >  a, 


(5.1) 


where  b  and  a  are  respectively  the  lengths  of  the  semi-major  and  minor  axes.  The  water 
has  density  pw  and  sound  speed  cw,  and  the  material  of  the  cavitator  is  assumed  to  be 
locally  reacting,  such  that  at  a  point  on  Sn  where  the  acoustic  pressure  is  p  and  the 
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outward  (into  the  water)  normal  component  of  velocity  is  vn  we  can  write 

—  Pw^w^  1  (5-2) 

Vn 

where  Z  is  the  dimensionless  surface  impedance.  In  practice  the  Mach  number  of  the 
relative  mean  flow  over  the  supercavity  is  sufficiently  small  that  the  convection  of  sound 
can  be  ignored.  Also  ignoring  turbulence  noise  sources  in  the  mean  flow,  the  acoustic 
pressure  p(x,  uj)elut  of  radian  frequency  uj  in  the  water  propagates  like  [6,  10,  12,  18,  32] 

(V2  +  k2w)p  =  0,  (5.3) 


where  kw  =  uj/cw  is  the  acoustic  wavenumber.  The  visco-thermal  dissipation  of  the  sound 
is  ignored,  so  that  equation  (5.2)  implies 


dp  ikwp 
dxn  Z 


0  on  SN, 


(5.4) 


where  xn  is  the  local  outward  normal  coordinate. 

Consider  the  sound  generated  in  the  water  by  the  impingement  of  ventilating  gas 
jets  on  the  cavity  wall  Sc-  The  jets  produce  a  localized  fluctuating  surface  pressure 
pv(x,ui)  on  Sc  which  is  equivalent  to  an  array  of  dipole  acoustic  sources  radiating  into 
the  water.  The  distribution  and  frequency  of  pv(x,uj)  are  known  from  the  measurements 
described  in  Chapter  4.  The  disparity  between  the  mean  gas  and  water  densities  allows 
the  interface  to  be  treated  as  a  pressure  release  surface,  and  we  therefore  assume  that 
pressure  fluctuations  elsewhere  on  Sc  can  be  neglected. 

Let  Pn{uj)  denote  the  acoustic  pressure  produced  by  jet  impingement  at  the  tip  of 
the  nose  Sn  (the  intersection  of  the  negative  x  axis  and  Sn).  To  express  this  in  terms  of 
pv(x,u)  we  introduce  a  Green’s  function  G(x,uj)  with  outgoing  wave  behavior  which  is 
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defined  by 

(V  2  +  kl )  G  =  0 

in  the  water,  subject  to 

+  =  Kv)Kz)  on  SN, 

G  =  0  on  Sc, 


(5.5) 


(5.6) 


where  S(-)  is  the  Dirac-delta  function  and  the  coordinates  y,z  are  defined  such  that 
x  =  (x,  y,  z )  forms  a  right-handed  Cartesian  basis.  This  defines  G  to  be  an  axisymmetric 
function  G(x,r,u>)  of  x  and  r  =  \Jy2  +  z2. 

Applying  Green’s  theorem  and  the  radiation  condition  (see  Appendix  B)  to  (5.5)  and 
(5.3)  yields  [2,  10,  18] 


/ 

dsN+sc 


p(x,  O')  W— -(x,  u)  —  G(x,  u>)J^-(x,  u) 


dxr 


dxr 


dS 


0, 


(5.7) 


where  the  integration  is  over  the  combined  surfaces  Sn  +  Sc  of  the  cavitator  and  the 
cavity  interface.  The  boundary  conditions  (5.4)  and  (5.6),  and  the  requirement  that 
p  =  pv  on  Sc  then  imply  that 


PnM 


p„(x,  u)  T(x,  u)  dS, 


(5.8) 


where  the  integration  is  over  the  impingement  region  S,  say,  of  the  cavity  interface  Sc, 
and  T(x,u)  is  a  transfer  function  defined  by 

T(x,w)  =  -  .  (5.9) 

\  /  r—a 

T(x,u)  depends  only  on  frequency  u>  and  distance  from  the  cavitator  edge  x;  the  ax- 


65 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


isymmetry  of  the  cavitator  and  cavity  demand  a  solution  for  the  acoustic  pressure  at 
the  nose  which  is  independent  of  the  azimuthal  location  of  the  point  on  Sc  where  pv  is 
applied. 


5.2  Calculation  of  the  transfer  function 


A 


Figure  5.2.  Reciprocal  problem  of  the  solution  of  the  wave  equation  for  the  motion 
produced  in  the  water  by  a  unit  point  source  at  the  nose  used  to  calculate  the  transfer 
function 

To  calculate  the  transfer  function  we  consider  the  reciprocal  problem  in  which  we 
first  calculate  the  motion  produced  on  the  jet  impingement  region  S  from  a  point  source 
at  the  nose  (Figure  5.2).  The  region  S  lies  in  the  geometric  shadow  of  direct  acoustic 
radiation  from  the  source.  The  behavior  of  the  Green’s  function  in  this  region  is  therefore 
governed  by  the  diffraction  of  sound  waves  around  the  nose  of  the  cavitator,  Sn  (see 
[8,  15,  24,  25,  31,  35,  37,  38,  47]).  This  will  be  analyzed  by  interpolation  between 
analytical  predictions  of  G(x,  r,  a;)  for  large  and  small  values  of  kwa ,  an  extension  of  the 
method  of  Howe  et  al.  [23] . 
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5.2.1  Surface  diffraction  when  kwa  1 

The  surface  diffraction  of  high  frequency  waves  from  the  source  can  be  represented  in 
terms  of  creeping  modes.  The  incidence  of  these  modes  on  the  edge  A  of  the  cavitator  is 
first  calculated;  a  separate  calculation  is  then  performed  to  determine  their  subsequent 
propagation  over  the  pressure  release  surface  Sc  to  impingement  region  S. 

The  creeping  wave  representation  is  found  by  choosing  the  local  approximation  G  = 
Gn,  where  in  the  immediate  vicinity  of  N  Gn  satisfies  the  ’earth  flattening’  approximation 
of  the  wave  equation  [8,  15,  25,  31,  35,  37] 


xn  <C  a,  (5.10) 


subject  to 


(5.11) 


where  w  denotes  curvilinear  ‘radial’  distance  on  Sn  from  the  point  source  at  the  nose 
along  a  meridional  arc  of  the  ellipsoid  (5.1),  and  p  =  a2 /b  is  the  radius  of  curvature  of 
Sn  at  N. 


Following  the  method  of  Howe  et  al.  [23],  we  introduce  the  zeroth-order  Hankel  trans¬ 
form 


(5.12) 


where  J0  is  the  Bessel  function  of  zero  order.  Therefore  [42] 


(5.13) 


where,  from  (5.11), 


dxn  Z  (27t)2’ 


(5.14) 
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The  solution  for  this  system  with  outgoing  wave  behavior  is  [8,  37] 


GN(xn,k)  = 


(2 


7T 


Ai 


x 


( m2  v 


Xr 


+  (§-f 


J(k-a 7-f) 


(5.15) 


m2 

P 


■Ai' 


,  L_f 


%Ai  <1  e~f  (  3  ~  1 


where  Ai  denotes  the  Airy  function  [1]  and  the  prime  denotes  differentiation  with  respect 
to  its  argument.  Putting  xn  =  0  on  the  surface,  we  find  from  (5.12)  that 


Gn(0,  w) 


1 

27T 


(5.16) 


x 


kJo(kzu)  Ai  <j 

[e-f  (^)§ 

(S-i)} 

ei(fcw-f)  dk 

[^(^(tM 

|  +  ^Ai 

{»-*  (^)f(S-l)} 

When  kwa  1  the  Bessel  function  may  be  replaced  by  its  large  argument  approximation 
( J0(x )  ~  a/2/ttx  cos(x— 7t/4)  [1]),  and  the  range  of  integration  extended  to  — oo  <  k  <  oo 
to  give 


Gn(0,  w) 


(2vr 


3 

12-x/tU 


(5.17) 


r»OQ 


yjk  +  ie  Ai 


x 


J(kru-f)  dk 


-oo  «^e-fAi/{e-f  (Me')1  (%-i)  j>  +  %iAi<|e-T  (  §  -  1 


P 


k2 


where  \k\  =  \Jk2  +  e2  (e  — >  +0),  with  the  branch  cuts  for  y/k  ±  ie  taken  respectively 
from  =p ie  to  =(=foo  in  the  /c-plane. 

As  kww  — >  oo  the  main  contribution  to  the  integral  is  from  the  residues  of  poles  at 
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the  zeros  k  =  kn  (n  =  1,  2,  3,  . . .)  of 


m 

p 


12L  A  ./ 

3  Ai 


tkw 


Ai 


/  \k\p 


=  0 


(5.18) 


in  the  upper  /c-plane.  By  making  use  of  the  formula  Aiw(£)  =  £Ai(£)  [1]  we  accordingly 
find 


Gn(0,  m) 


‘Mjt 

3e  j 


E 


\fknZ 


ik-n 


el\ph TW  f  kw  _  iknp 

"  Zkn 


iknp 

2 


2k  w  I  1 
L.2  I  -L 


i-ii-i)! 


kwa  1.  (5.19) 


The  solution  is  valid  for  kQw  1  and  is  asymptotically  accurate  for  w  less  than  about 
7 a  [35].  Howe  et  al.  [23]  demonstrate  how  the  geometrical  theory  of  diffraction  [31,  37] 
can  be  used  to  extend  the  range  of  validity  of  the  solution  out  to  the  edge  A  of  the 
cavitator.  To  do  this  the  algebraic  amplitude  factor  \j\fm  involving  the  curvilinear 
surface  propagation  distance  m  must  be  replaced  by  1  /y/r  and  the  complex  phase  knw  in 
the  exponential  factor  elknm  must  be  modified  to  account  for  the  variation  of  the  creeping 
wavenumber  kn  =  kn{m )  because  of  the  change  in  the  radius  of  curvature  p  =  p(m)  of 
creeping  surface  rays  with  distance  m  along  the  propagation  path  on  Sn-  The  result  can 
be  expressed  in  the  form 


Gn(0,  m) 


3i7T 

3e  4 


E 


\fkn 


2v/2vr r  ^  \  Al_ 

"  '  Zkn 


o  I  1.2  T  1 

""n 


(i-*: 


where  the  quantity  (-)o  is  evaluated  at  the  nose  N  (where  p  = 
given  by 


-r  e1^,  kwa  »  1. 

Si  / 

kl  \  /0 

(5.20) 

a2 /b)  and  the  phase  is 


r>zn 

=  /  kn(m’)  dm', 
Jo 


(5.21) 
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in  which  kn  is  the  appropriate  solution  of  Equation  (5.18)  when  p  =  p{w'). 

5.2.2  Diffraction  at  the  cavitator  edge 

When  the  creeping  surface  wave  (5.20)  reaches  the  edge  A  of  the  cavitator  (Figure  5.2) 
it  is  further  diffracted  because  of  the  change  in  surface  condition  from  (5.11)  on  Sn  to 
G  —  0  on  Sc-  Near  A  (where  r  — »  a)  the  approaching  surface  wave  Gn(0,  w)  has  a  slowly 
varying  amplitude  but  rapidly  varying  phase  factor  ~  e*fc™TO.  The  diffraction  of  this  wave 
at  the  edge  was  treated  approximately  by  Howe  et  al.  [23]  by  imagining  the  surface  of 
the  cavitator  to  be  deformed  near  A  into  a  circular  cylinder  of  radius  a  (depicted  by 
the  broken  lines  in  Figure  5.2)  that  forms  a  rigid  upstream  continuation  of  the  mean 
cylindrical  cavity  interface  Sc-  The  creeping  wave  (5.20)  is  then  regarded  as  incident  on 
A  from  along  this  cylinder,  and  expressed  in  the  form 

Gn(0,  w)  —  QieikwX  (5.22) 


where 


Qi  = 


o  Sit 

3e  4 


£ 


2v/2vra  ^  \  Jsml. 


VK, 


ikna 2 
2  b 


I  1 

7.2  I  T 


^47) 


(5.23) 


and  £  =  bE  b ^  J  curvilinear  distance  on  the  meridional  geodesic  on  the 

ellipsoidal  cavitator  Sn  from  the  nose  N  to  the  edge  A,  where  E(£)  is  the  complete 
elliptic  integral  of  the  second  kind  [1], 

Howe  et  al.  [23]  solved  this  diffraction  problem  for  an  acoustically  hard  cavitator  by 
the  method  of  Weiner-Hopf  [9,  33]  in  the  limiting  case  where  kwa  1;  the  same  method 
is  applicable  here.  Letting  Gg  denote  the  edge  diffracted  component  of  G ,  in  the  region 
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r  >  a  we  can  write 


G  —  QieikwX  +  Gg(x,  r), 


(5.24) 


where  for  some  function  A(k)  we  can  put  [10,  18,  33] 


/  LXJ 

Gs(x,r)=  A(k)H{01]Meikxdk, 


(5.25) 


where  Hq1'1  is  the  zeroth  order  Hankel  function  and  7  =  \Jk\  —  k2,  with  branch  cuts 
in  the  /c-plane  extending  from  k  =  ±(kw  +  ie)  (e  — >  +0)  to  ±00,  so  that  7  is  positive 
imaginary  on  the  real  axis  when  \k\  >  kw.  The  function  A(k)  is  determined  from  the 
conditions 


dG  I  ikwG  _  Q  r  ^  0 
W  +  ~Z~  ~  Uj  x<[) 


,  r  =  a  +  0 


(5.26) 


G  =  0,  x  >  0 

which  supply  the  following  dual  integral  equations  satisfied  by  A{k) 


'yH[1\iya)  +  —^-H^\^a)  ]  A(k)elkxdk  =  0,  x  <  0,  (5-27) 


Z 


A(k)H^(la)  + 


2m  (k  —  kw  —  ie) 


eikxd k  =  0,  x  >  0. 


(5.28) 


For  kwa  3>  1,  the  ffankel  functions  Hq  ,  can  be  replaced  by  their  large  argument 
approximations  [1],  Solving  by  the  Wiener- ffopf  method  [9,  33] 


i  7  r  K+(kw)e‘^+^-“»dk 


K+(k){k  —  kw  —  ie) 


,  kwa  1  (5.29) 


where  K+(k)  is  one  of  two  functions  K±(k)  regular  on  the  real  k  axis  and  respectively 
in  Im  k  <  0,  and  satisfying  for  real  k 


k. 


K+(k)K-(k)  =  -G-  +  \J (kw  +  ie)2  -  k 2 


Z 


(5.30) 
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Taking  the  derivative  of  (5.29)  on  the  cavity  interface  r  =  a, 


f <9G\  _  Qi  f°°  K+(kw)\/kw  +  ie  +  k 

\dr  ) r=a  2vr  J_x  K+(k)y/kw  +  ie  -  k 


e  — »■  +0,  kwa  1, 


(5.31) 


When  x  >  0  and  kwx  is  large,  the  main  contribution  to  the  integral  in  (5.31)  is  from  the 
neighborhood  of  the  branch  point  at  k  =  kw  +  ie.  By  expansion  of  the  integrand  about 
this  point  we  find,  to  leading  order,  that 


dG 

dr 


-tH  ~  Qi 


2  k 


nx 


—  ei(kwX  4)  kwa^>  1. 


(5.32) 


Using  (5.23)  and  the  definition  (5.9),  we  obtain  from  this  the  corresponding  explicit 
approximation  for  the  transfer  function 


T(x,u)  ~ 


\/  kna 


ikna2  f  1 

2  b  [  kl  +  1 


Jl'k  n(£)  +  kwx} 


(5.33) 


for  kwa  1. 


5.2.3  Behavior  of  the  transfer  function  when  kwa  ~  0(1) 

The  approximate  behavior  of  the  transfer  function  when  kwa  ~  0(1)  is  obtained  by 
the  method  discussed  by  Howe  et  al.  [23]  in  which  the  surface  wave  on  Sn  incident  on 
the  edge  A  of  the  cavitator  is  approximated  by  the  surface  potential  produced  by  the 
point  source  at  N  on  a  flat ,  infinitely  extensive  surface.  In  a  first  approximation  at  low 
frequencies,  this  potential  is  the  same  as  that  generated  by  a  source  on  a  rigid  surface, 
for  which 

zu 

Gn(0,w)  =  -- - .  (5.34) 

Znzu 
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The  influence  on  this  wave  of  diffraction  at  the  edge  A  will  be  assumed  to  be  well 
approximated  by  the  formula  (5.32)  of  §5.2.2,  wherein  we  take 


r±ikwd. 

e'  =  -  2 S* 


(5.35) 


where  £  is  defined  as  in  equation  (5.23).  Therefore  (5.32)  yields 


T{x,u)~  t  3  kwa~0{  1).  (5.36) 

22iT2a£  y  x/a 


5.3  Transfer  function:  Numerical  results 


The  properties  of  a  compliant  cavitator  are  most  easily  demonstrated  by  letting  Z  —  1 
in  the  surface  condition  (5.4)  on  Sn,  corresponding  to  a  locally  reacting  Lpc  surface’.  The 
attenuation  of  creeping  waves  on  Sn  is  then  governed  by  the  solutions  kn  of  equation 
(5.18))  when  Z  —  1.  At  very  high  frequencies  {kwp  — >  oo)  the  zeros  are  the  same  as 
those  for  a  pressure  release  surface  [37],  which  are  determined  by  the  roots  £  =  —  £n,  n  = 
1,  2,  3,  ...  of  Ai(£)  =  0,  which  are  real  and  negative  and  given  by  [1] 

2 

3  ,  n  >  3.  (5.37) 


fi  =  2.33818,  &  =  4.08794,  £3  =  5.52056,  = 


—  (in  — 

8  V 


Each  £n  determines  a  corresponding  pole  kn  given  approximately  by 

kn  £n  (  2  \  ®  in  .  . 

-T  ~  1  +  T  7—  e^’  n  =  1,  2,  3,  ■  ■  ■ ,  kwp  »  1.  (5.38) 

kw  2  y  kw  p  J 

The  exact  dependencies  on  kwp  of  the  real  and  imaginary  parts  of  the  solutions  kn/kw 
of  equation  (5.18)  are  plotted  in  Figures  5.3,  5.4  for  n  —  1,  2,  3.  Re(fcn)  ~  kw  when  kwp 
exceeds  about  10,  and  Im(fcn)  increases  rapidly  with  n  for  fixed  kwp.  Indeed,  as  a  result, 
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it  is  readily  verified  that  only  the  first  term  (n  —  1)  makes  a  significant  contribution  to 
the  sum  on  the  right  hand  side  of  (5.33). 


l  to  100  1000 

k  w 

Figure  5.3.  Real  parts  of  the  creeping  mode  wavenumbers  kn/kw  for  n  1.  2,  3  for  a 
pc- cavitator  (Z  =  1) 

The  comparison  in  Figure  5.5  of  the  numerical  solution  of  (5.18)  and  the  asymptotic 
approximation  (5.38)  for  this  principal  creeping  mode  (n  =  1)  indicates  that  (5.38)  may 
be  used  when  kwp  >  10.  In  particular  it  can  be  used  to  evaluate  the  phase  T1(£)  for  the 
first  term  in  the  expansion  (5.33).  To  do  this  it  is  convenient  to  introduce  the  parametric 
representation  (x,  r)  =  (— b cos  (ft,  asin0),  0  <  0  <  of  a  point  on  the  ellipsoidal  nose 
(5.1).  Then 


P  = 


L  /  2  2 
—  (a  cos 

ab 


0  +  b2  sin2  0) 


dw  = 


(2  2 
(a  cos 


0  +  b2  sin2  0)  2  d0, 


(5.39) 
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Figure  5.4.  Imaginary  parts  of  the  creeping  mode  wavenumbers  kn/kw  for  n  =  1,  2,  3 
for  a  pc- cavitator  (Z  =  1) 


and 


*i( 


~  kw£  + 
=  kwl  + 


kwa 

~2b 

knna 


e  3 


dcj) 


'o  \/ sin2  (j)  +  ( a2/b 2)  cos2  <f> 


2  X  l 


2b 


K 


b2  -  a 2 
b2 


e  3  , 


(5.40) 


where  K  denotes  the  complete  elliptic  integral  of  the  first  kind  [1]. 

The  magnitude  of  the  transfer  function  \T(x,uS)\  determines  the  strength  of  the  self 
noise  at  the  nose  N  in  terms  of  the  impinging  gas  pressure  held  pv  on  the  cavity  interface. 
Its  variation  for  kwa  1  is  given  by  the  first  mode  in  (5.33)  in  the  form 


a2\T(x,cu)\ 


3 

2tt 


( k±a )2e*'I'1W 


kw  _  ik\a-  f  2fcm  .  -i 

fcl  2 6  \  fcf  ~l~ 


(5.41) 


This  approximation  is  plotted  as  the  ‘pc-nose’  curve  in  Figures  5. 7-5. 9  for  a  cavitator 
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Figure  5.5.  Comparison  of  the  real  and  imaginary  parts  of  k\/kw  for  a  pc  cavitator 
determined  by  equation  (5.18)  (with  Z  =  1)  and  the  asymptotic  approximation  (5.38) 
with  ri—1. 

of  aspect  ratio  b/a  =  5/3  and  x  =  a,  2a,  3a  downstream  of  the  edge  A  of  the  cavitator. 
The  straight  line  labeled  ‘flat  nose’  in  the  figure  is  the  corresponding  low  frequency 
prediction  given  by  (5.36)  (in  which  l/a  =  (b/a)E( 0.8)  ~  2.125).  The  behavior  at 
intermediate  frequencies  is  estimated  by  a  smooth  interpolating  curve  (—  —  — )  between 
these  high  and  low  frequency  limits. 

Figures  5. 7-5.9  also  depict  the  corresponding  prediction  of  | T(x,  u ;)  |  for  an  acoustically 
hard  ellipsoidal  cavitator  of  the  same  aspect  ratio  and  distances  from  the  cavitator  edge, 
which  obtains  in  the  formal  limit  Z  — >  oo.  In  this  case  equation  (5.18)  for  the  creeping 
mode  poles  reduces  to 

Ai'{e-¥  (ihC)1  (§-i)}=°.  (,42) 

Again  only  the  pole  k\  in  the  upper  half  plane  nearest  to  the  real  k  axis  makes  a  significant 
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Figure  5.6.  Comparison  of  the  real  and  imaginary  parts  of  k\jkw  for  an  acoustically 
hard  cavitator  ( Z  — >  oo)  with  the  approximation  (5.43)  with  n  =  1. 


contribution  to  the  creeping  wave  amplitude.  The  calculated  dependence  of  the  real  and 
imaginary  parts  of  k\/kw  on  kwp  are  shown  in  Figure  5.6.  The  asymptotic  approximation 
is  given  by 

k-t  STi  f  9  \  3  in 

eT,  cr  i  =  1.01879.  (5.43) 


h  ,  ,  (J\ 

—  ~  1  4 - 

h  9 


kwp 


The  zeros  of  Ai^t))  lie  on  the  negative  real  axis,  and  (  =  —  Ui  is  the  root  closest  to  the 
origin.  The  approximation  (5.43)  has  been  used  to  evaluate  the  complex  phase  T1(£)  for 
the  hard  cavitator,  yielding 


'IT  (A)  =  kw£  +  <Ti 


(5.44) 
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Figure  5.7.  Frequency  dependence  of  the  creeping  wave  and  low  frequency  (‘flat  nose’) 
approximations  for  20  x  log10(a2|T|)  (dB)  at  x  =  a  for  an  ellipsoidal  cavitator  of  aspect 
ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc- nose  or  is  acoustically  hard;  the 
broken  line  curves  represent  possible  interpolations  for  intermediate  frequencies. 

5.4  Prediction  of  the  sound  pressure  spectrum  at  the 
nose 


To  predict  the  self-noise  at  the  nose  of  our  experimental  supercavitating  vehicle  (Fig.  1.1) 
we  consider  that  the  fluctuating  surface  pressure  pv  is  generated  by  an  axisynnnetric  array 
of  20  jets  of  infinitesimal  cross-section  discreetly  impinging  the  surface  Sc,  hence 

Pv  =  y^Fn(x,  u)5(x  -  xn)— — —  (5.45) 

z ^  n 
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Figure  5.8.  Frequency  dependence  of  the  creeping  wave  and  low  frequency  (‘flat  nose’) 
approximations  for  20  x  log10(a2|T|)  (dB)  at  x  =  2a  for  an  ellipsoidal  cavitator  of  aspect 
ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc- nose  or  is  acoustically  hard;  the 
broken  line  curves  represent  possible  interpolations  for  intermediate  frequencies. 


where  Fn  is  the  unsteady  force  exerted  by  the  nth  jet  applied  at  its  centroid  ( xn ,  6n )  and 
6  is  the  polar  angle  about  the  x-axis.  Using  this  to  approximate  the  integral  (5.8), 

20 

Pn(uj)  =  Y 

n=  1 
20 

=  YFn(xn,On,uj)T(xn,uj).  (5.46) 

n= 1 

where  T(xn,  tv)  is  the  transfer  function  corresponding  with  the  centroid  of  the  nth  jet.  As 
we  are  interested  in  the  spectral  content  of  the  noise  at  the  nose  &p,n  =  ^{Pn{u)p*n{uj)), 
we  multiply  each  side  by  its  complex  conjugate  and  divide  by  2i r  to  find  an  expression 


Fn(x,  uj)5(x  —  xn)5{6  —  9n)T(x,u>)  dxdd 
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Figure  5.9.  Frequency  dependence  of  the  creeping  wave  and  low  frequency  (‘flat  nose’) 
approximations  for  20  x  log10(a2|T|)  (dB)  at  x  =  3a  for  an  ellipsoidal  cavitator  of  aspect 
ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc- nose  or  is  acoustically  hard;  the 
broken  line  curves  represent  possible  interpolations  for  intermediate  frequencies. 


for  the  sound  pressure  spectrum  at  the  nose 
1  20  20 

$p,jvM  =  —  '^2'^2(Fn(xn,6n,u})F^l(xn,6n,u)))T(xn,u))T*(xm,u))  (5.47) 

71=1  771=1 

where  *  denotes  the  complex  conjugate.  As  the  jets  are  assumed  to  be  statistically 
independent,  the  terms  of  the  double  summation  are  non-zero  only  if  n  —  m,  leaving 

20 

$P,n(v)  —  $ FF,n(w)\T (xn,  u)\2  (5.48) 

71=1 

where  $>FF,n{ w)  is  the  frequency  spectrum  of  the  force  exerted  by  the  nth  jet  on  the 
surface.  If  we  consider,  as  in  Chapter  3,  that  the  20  jets  are  statistically  independent 
and  that  they  are  equidistant  from  the  nose  so  that  X\  =  X2  =  X3...X20  =  x,  then  (5.48) 
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reduces  to 

®p,n(w)  =  (w)  |  T(x,uj)  | 2  (5.49) 


5.5  Spectrum  predictions  at  the  nose 

The  experimental  results  of  §4.4  and  the  numerical  results  of  §5.3  are  now  used  in 
Eq.  (5.49)  to  predict  the  impingement-generated  sound  pressure  spectrum  at  the  nose 
of  the  experimental  supercavitating  vehicle.  Figures  5.10-5.12  present  the  normalized 
prediction  for  the  sound  pressure  spectra 


10  x  log10 


(  ®P,n(w)  \ 
\Mp2AU3a  ) 


versus  dimensionless  frequency  kwa  for  jets  located  at  distances  x/a  =  1,  2,  and  3  aft 
of  the  cavitator  edge,  where  M  —  U/cw  is  the  Mach  number  of  the  jet  relative  to  the 
sound  speed  in  water  cw,  Pa  is  the  mean  jet  density  (taken  as  1.2  kg/m3)  and  U  is  the 
jet  impact  velocity.  Predictions  for  both  the  hard  and  compliant  (pc)  nose  are  given 
using  the  experimentally  measured  for  a  gas  injection  rate  of  Q  =  0.0094  m3/s 

and  i  =  35  mm,  chosen  for  being  the  most  typical  parameters  of  ARL’s  model  vehicle 
experiments.  U  is  therefore  taken  as  31.1  m/s  (see  Table  4.2).  Predictions  made  using  the 
force  spectrum  measurements  from  Chapter  4  for  each  measured  jet  length  and  injection 
flow  rate  are  presented  in  Appendix  E. 

Differences  between  sound  at  the  nose  for  the  hard  and  pc  nose  are  small  at  low 
frequencies  [kwa  <  1),  but  very  large  at  high  frequencies  ( kwa  ~  10),  on  the  order  of 
30dB  or  more.  The  sound  spectra  all  peak  near  kwa  =  0.1,  with  the  predictions  for 
a  hard  nose  dropping  off  like  (/c^a)-1  and  the  predictions  for  the  pc  nose  dropping  off 
like  ( kwa)~ 2 .  There  is  approximately  a  2dB  decrease  in  the  sound  amplitude  between 
the  predictions  when  the  jets  are  one  and  two  cavity  radii  away  from  the  cavitator  edge 
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Figure  5.10.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a  for  an  ellipsoidal 
cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose  or  is 
acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length  35  mm. 


(x  =  a,  2a),  and  a  1.5dB  decrease  between  predictions  for  jets  two  and  three  cavity  radii 
away  (x  =  2a,  3a). 


5.6  Conclusions 

Comparing  the  predictions  for  the  sound  at  the  nose  with  the  predictions  for  the  radiated 
sound  in  the  water  found  in  Chapter  4,  there  is  a  dramatic  decrease  in  the  presence  of 
high  frequency  sound  at  the  nose,  particularly  for  a  compliant  (pc)  nose  section.  This 
suggests  that  the  high  frequencies  (kwa  >  10)  generated  by  the  jet  impingement  source 
have  very  low  amplitudes  at  the  nose  and  will  have  negligible  impact  on  guidance  systems 
located  forward  of  the  cavitator. 

Although  we  predict  that  high  frequencies  at  the  nose  will  be  minimal,  it  is  always 
desirable  to  reduce  the  self-noise  generated  in  and  around  the  supercavity.  An  important 
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Figure  5.11.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  2 a  for  an  ellipsoidal 
cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose  or  is 
acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length  35  mm. 


conclusion  from  our  predictions  for  the  sound  spectrum  at  the  nose  which  applies  to  the 
self-noise  problem  is  that  high  frequencies  are  significantly  damped  by  an  acoustically 
soft  nose;  we  predict,  for  example,  a  30  dB  decrease  in  high  frequencies  between  the 
hard  and  compliant  nose  sections.  Furthermore,  small  decreases  to  the  amplitude  of  the 
broadband  sound  at  the  nose  may  be  achieved  by  increasing  the  distance  between  the 
edge  of  the  cavitator  and  the  ventilating  jets.  This  insight  provides  two  strategies  for 
reducing  the  high  frequency  noise  at  the  nose  of  a  supercavitating  vehicle.  The  first  is 
to  choose  the  most  acoustically  compliant  material  for  the  nose  section  which  still  allows 
for  transmission  of  guidance  communications  and  does  not  inhibit  the  vehicles  hydrody¬ 
namic  performance.  The  second  is  to  locate  the  ventilation  jets  as  far  downstream  from 
the  cavitator  as  possible  without  inhibiting  the  generation  of  the  supercavity  or  vehicle 
performance. 
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Figure  5.12.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  3a  for  an  ellipsoidal 
cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose  or  is 
acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length  35  mm. 
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Chapter  6 
Conclusions 

6.1  Summary 

Ventilated  supercavities  yield  dramatic  reductions  in  vehicle  drag,  but  they  can  be  a 
significant  source  of  noise,  some  of  which  can  inhibit  successful  deployment  of  a  super- 
cavitating  underwater  vehicle.  The  research  described  in  this  dissertation  was  motivated 
by  the  need  to  understand  the  contributions  to  this  noise  made  by  the  ventilating  jets  of 
the  supercavity.  The  problem  was  approached  by  the  development  of  three  general  theo¬ 
retical  models  for  jet-cavity  interaction  supported  by  experiments.  Empirical  knowledge 
of  the  cavity  and  ventilating  jets  determined  from  a  scale-model  experimental  vehicle 
were  used  with  the  theory  to  make  specific  predictions  for  the  sound  and  self-noise  gen¬ 
erated  by  the  ventilating  jets. 

The  Erst  theoretical  problem  involved  an  idealized,  gas-hlled  spherical  cavity  whose 
interface  was  excited  by  a  planar  ring  of  radially  projecting  jets  originating  from  a  source 
at  the  sphere’s  center.  Analysis  of  this  configuration  provided  a  general  understanding 
of  the  directivity  of  the  sound  radiated  into  the  water  as  well  as  the  importance  of 
cavity  resonances  and  damping  mechanisms  within  the  cavity,  ft  was  deduced  that 
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low  frequency  sound  emanates  from  the  spherical  cavity  omni-directionally,  while  higher 
frequencies  are  generated  by  a  source  of  dipole  type,  projecting  sound  outward  and 
normal  to  the  cavity  at  the  jet  impingement  site.  It  was  determined  that  the  amplitudes 
of  resonant  cavity  oscillations  are  controlled  by  sound  radiation  damping  and  by  thermal 
and  viscous  losses  within  the  cavity;  it  was  argued  that  the  steady  exhausting  of  gas 
at  the  rear  of  a  real  supercavity  has  a  negligible  effect  on  damping.  Another  important 
conclusion  was  confirmation  of  the  validity  of  the  effective  pressure  release  nature  of  the 
gas-water  interface,  which  implies  that  the  dominant  source  of  sound  radiated  into  the 
water  is  associated  with  the  direct  impact  of  the  ventilating  jets  on  the  cavity  interface. 

The  second  theoretical  problem  examined  the  sound  generated  in  the  water  by  a 
turbulent  gas  jet  of  infinitesimal  cross-section  impinging  at  normal  incidence  on  a  plane 
gas- water  interface.  The  acoustic  spectrum  was  represented  in  terms  of  the  unsteady  force 
exerted  on  the  interface  by  the  jet.  Theoretical  results  were  combined  with  measurements 
of  the  force  spectrum  for  a  jet  impinging  on  a  rigid  wall  to  predict  the  acoustic  spectrum 
in  the  water.  Results  for  a  plane  interface  were  then  applied  to  develop  an  understanding 
of  the  sound  generated  for  a  more  realistic,  cylindrical  cavity. 

Finally,  a  theory  was  developed  to  estimate  the  high-frequency  self-noise  at  the  vehicle 
nose  in  terms  of  the  pressure  fluctuations  at  the  cavity  interface  induced  by  the  ventilating 
jets.  This  was  done  by  deriving  an  asymptotic  approximation  for  a  transfer  function  that 
accounts  for  the  diffraction  of  sound  around  an  ellipsoidal  surface  of  the  vehicle  nose. 
Numerical  predictions  of  the  transfer  function  revealed  the  important  role  played  by  the 
surface  impedance  of  the  nose;  an  acoustically  hard  nose  results  in  significantly  greater 
acoustic  amplitudes  at  the  nose  than  an  acoustically  soft  or  compliant  nose. 

It  was  determined  that  as  a  first  approximation  the  supercavity  interface  could  be 
treated  as  rigid.  This  allowed  for  empirical  knowledge  of  the  ventilating  jets  of  an 
experimental  supercavitating  vehicle  to  be  acquired  by  means  of  an  experiment  involving 
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the  measurement  of  the  frequency  spectrum  of  the  unsteady  force  exerted  by  a  gas  jet  on 
a  plane,  rigid  surface.  Results  were  combined  with  the  theoretical  predictions  to  make 
specific  predictions  of  the  frequency  dependence  and  amplitude  of  the  sound  generated  by 
a  supercavitating  vehicle.  The  experiments  also  confirmed  an  important  approximation 
of  the  theoretical  model,  namely  that  the  impingement  source  region  is  acoustically 
compact.  The  predictions  of  the  sound  generated  by  the  measured  supercavity  jets 
were  compared  with  predictions  made  using  an  independent  measurement  of  the  surface 
pressure  on  a  flat  plate  produced  by  an  air  jet  impinging  at  normal  incidence.  The 
close  coincidence  between  these  predictions  lends  valuable  support  for  the  validity  of  our 
measurement  procedure. 

6.2  Recommendations  for  future  work 

6.2.1  Extended  experiments 

The  work  described  in  this  dissertation  compliments  experimental  investigations  of  the 
acoustics  of  ventilated  supercavities  on  underwater  vehicles.  Our  predictions  for  the  self 
noise  at  the  nose  should  be  investigated  by  the  direct  measurement  of  the  sound  at  the 
nose  of  a  deployed  supercavitating  vehicle  that  uses  the  modified  gas  injection  system 
with  radially  projecting  ventilation  jets.  Our  results  and  the  advanced  acoustic  levels 
produced  by  the  impact  of  the  jets  on  the  cavity  interface  at  normal  incidence  will  help 
isolate  impingement-generated  sound  from  that  created  by  other  acoustic  sources  within 
and  around  the  cavity.  The  experiment  should  be  conducted  using  a  free-running  vehicle 
rather  than  one  mounted  in  a  water  tunnel,  where  the  enclosed  test  section  creates  a 
reverberant  field  that  may  influence  the  measured  self-noise  levels. 

As  the  experimental  supercavitating  vehicles  evolve  in  scale  and  design,  it  will  be 
desirable  to  make  predictions  for  the  impingement-generated  sound  specific  to  each  velii- 
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cle.  The  theoretical  solutions  are  general  and  applicable  to  any  cavity  with  a  comparable 
injection  system.  If  acoustic  compactness  of  the  impingement  source  is  maintained,  that 
is,  if  the  jet  impact  diameter  remains  small  relative  to  the  wavelength  of  the  relevant 
acoustic  frequencies,  the  sound  predictions  are  easily  made  using  empirical  jet  and  cavity 
information  unique  to  each  vehicle.  The  benchtop  test  described  in  Chapter  4  can  be 
easily  setup  and  quickly  executed  using  either  a  vehicle  jet  or  a  stand  alone  jet  nozzle, 
facilitating  the  production  of  theoretical  results  which  compliment  the  acoustic  investi¬ 
gations  for  each  model  vehicle. 

A  possible  extension  to  the  self-noise  predictions  would  be  the  development  of  so¬ 
lutions  for  the  self-noise  generated  by  jets  impinging  at  shallow  incidence.  This  ar¬ 
rangement  more  closely  represents  the  acoustic  source  present  due  to  the  traditional  gas 
injection  system  on  supercavitating  underwater  vehicles  which  includes  gas  deflectors 
which  redirect  the  ventilating  jets.  The  prediction  method  could  again  be  derived  in 
terms  of  the  unsteady  force  exerted  by  a  jet  on  the  interface,  and  an  analogous  exper¬ 
iment  for  measuring  this  force  could  be  developed.  The  importance  of  considering  this 
more  general  problem  will  become  apparent  as  understanding  of  all  of  the  various  acous¬ 
tic  sources  in  the  supercavity  improves.  It  is  unlikely  that  a  theory  that  uses  a  more 
realistic  geometry  of  the  gas- impingement  source  would  reveal  any  new  insight  relevant 
to  the  problem  of  noise  reduction;  it  would  primarily  provide  a  more  accurate  prediction 
for  the  absolute  levels  of  the  impingement-generated  sound. 

6.2.2  Recommendations  for  achieving  noise  reduction 

This  dissertation  is  only  a  first  step  in  the  complete  investigation  of  the  self-noise  problem. 
However,  it  immediately  provides  valuable  insight  into  the  possible  means  of  reducing 
the  self-noise  at  the  nose  of  a  supercavitating  vehicle.  The  first  insight,  revealed  by  the 
measurements  of  Chapter  4,  is  the  relationship  between  ventilation  gas  injection  rate 
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and  the  amplitude  of  the  generated  sound.  A  doubling  of  the  ventilation  rate  produces 
a  rapid  increase  in  generated  sound  power,  approximately  10-15  dB.  This  leads  us  to 
recommend  an  investigation  to  determine  the  lowest  gas  injection  rate  capable  of  safely 
maintaining  a  stable  supercavity  around  the  underwater  vehicle. 

Chapter  5  demonstrated  that  the  impedance  of  the  cavitating  nose  section  material 
has  a  significant  impact  on  self-noise  at  the  vehicle  nose.  Predictions  for  an  acousti¬ 
cally  hard  cavitator  were  dramatically  louder  than  those  for  an  acoustically  compliant 
nose.  This  result  reinforces  the  importance  of  identifying  a  material  for  the  nose  which 
minimizes  the  diffraction  of  sound  from  its  trailing  edge  without  compromising  the  trans¬ 
mission  of  guidance  system  communications  or  hydrodynamic  performance. 

Another  important  conclusion  from  Chapter  5  is  the  demonstrated  effect  on  the 
sound  at  the  nose  of  the  distance  between  the  vehicle’s  ventilating  jets  and  the  edge 
of  the  cavitating  nose  section.  Proximity  of  the  jets  to  the  nose  significantly  increased 
the  magnitude  of  the  self-noise  at  the  nose.  We  therefore  recommend  exploring  the 
maximum  distance  between  the  vehicle’s  nose  section  and  gas  injection  site  for  which 
the  supercavity  is  successfully  and  safely  developed  and  maintained  without  adversely 
affecting  the  vehicle’s  hydrodynamic  performance  or  stability. 
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Appendix  A 

Review  of  select  topics  in  linear 
acoustics 


This  appendix  summarizes  some  of  the  fundamentals  of  acoustic  theory  pertinent  to 
this  dissertation.  The  information  herein  is  well  known  and  included  for  completeness. 
Material  is  taken  from  Blackstock  [5],  Howe  [18],  and  Pierce  [37].  See  these  references 
and  their  bibliographies  for  greater  detail  on  each  topic  or  a  more  general  review  of  the 
principles  of  linear  acoustics. 

A.l  The  linear  wave  equation 

The  acoustic  wave  equation  is  derived  from  the  conservation  and  constitutive  equations 
which  describe  a  fluid.  We  consider  a  fluid  volume  V  defined  by  a  surface  S  in  which 
body  forces  and  thermal  and  viscous  effects  are  negligible.  The  fluid  inside  the  volume 
has  velocity  u(x,  t),  pressure  P(x,  t),  and  density  p(x,  t)  at  a  position  x  =  (xi,x2,x3) 
and  time  t. 

Conservation  of  mass  requires  that  the  increase  of  mass  inside  V  must  be  equal  to 


90 


Report  No.  ME09  -  13 


Boston  University,  College  of  Engineering 


the  influx  of  mass  through  S.  This  is  expressed  by  the  continuity  equation 

^  +  u  =  0,  (A.l) 

where  the  material  derivative  D/Dt  =  d/dt  +  u  •  V. 

Conservation  of  momentum  holds  that  the  time  rate  of  change  of  momentum  inside 
V  is  equal  to  the  sum  of  all  forces  acting  on  the  fluid  volume  and  the  momentum  flux 
through  the  surface  S.  Neglecting  body  and  viscous  forces,  momentum  conservation  can 
be  expressed  by 

+  VP  =  0.  (A.2) 

Finally  we  require  an  equation  of  state,  a  relationship  between  thermodynamic  vari¬ 
ables.  We  will  assume  that  our  fluid  volume  is  isentropic,  that  is,  the  entropy  s  in  the 
volume  is  constant,  e.g. 

P  =  P(p)  (A.3) 

We  consider  that  the  pressure  P  and  density  p  can  be  decomposed  into  a  mean,  or  static, 
and  excess  constituent,  such  that 


:  Po  +  P, 

(A.4) 

Po  +  Sp, 

(A.5) 

where  pQ  and  pQ  are  respectively  the  mean  pressure  and  density  in  the  fluid  volume,  5p 
is  the  excess  density,  and  p  is  the  excess  or  acoustic  pressure.  This  leads  to  the  Taylor 
expansion  of  our  equation  of  state  (A.3) 


P 


Po  + A 


P~  Po 
Po 


B 

2! 


(W 


+  ... 


(A.6) 
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where  A,B ,  etc.  are  empirically  determined  coefficients.  We  then  introduce  the  sound 
speed  c  in  the  fluid,  defined  by 


c 


2 


dP 

dp 


A  +  B  P~  P°  + 

Po  Po  Po 


(A.7) 


Letting  p  — »  pa,  we  find  that  c  becomes  a  constant  cQ  and  that  A  =  p0c20 .  Applying  the 
decomposition  (A. 4,  A. 5),  we  arrive  at  the  equation  of  state 


p  =  cl  5  p 


(A.8) 


To  simplify  these  equations,  we  consider  that  fluctuations  in  the  fluid  properties 
induced  by  sound  energy  are  very  small.  In  a  quiet  fluid,  P  =  pQ,  p  =  pa,  and  u  =  0 
across  V.  We  also  consider  that  the  fluid  is  homogeneous,  and  that  pressure  and  density 
are  constant  across  the  volume.  We  will  assume  that  when  sound  is  introduced  the  excess 
pressure  and  density  and  induced  particle  velocity  are  very  small  compared  to  their  static 
values,  and  take  the  relations 


\Sp\  <  Po, 

(A.9) 

\p\  <  PoC20 , 

(A.  10) 

\u\  <C  ca. 

(A.  11) 

Applying  the  decomposition  (A. 5),  we  expand  the  continuity  equation  (A.l)  as 

d 

—  (8p)  +  uVSp  +  pQV  ■  u  +  <5pV  ■  u  =  0.  (A. 12) 

Considering  our  small  value  approximations,  we  note  that  the  second  and  fourth  terms  on 
the  left-hand  side  (LHS)  above  are  an  order  of  magnitude  smaller  than  the  first  and  third, 
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and  are  therefore  considered  negligible.  Eliminating  them,  we  arrived  at  the  linearized 
continuity  equation 

^(5p)+p0V-  u  =  0.  (A. 13) 

We  next  apply  the  same  procedure  to  the  conservation  of  momentum  equation  (A. 2). 
Expanding  (A. 2)  in  terms  of  the  decompositions  of  pressure  and  density, 


<9u  _  .  <9u  _ 

po yr  +  Ad*  ■  Vu  +  dp—  +  bpu  ■  Vu 
at  at 


Vp  =  0. 


(A.  14) 


As  before  we  consider  the  order  of  each  term,  and  find  that  only  the  first  and  fifth  terms 
on  the  LHS  are  of  first  order.  Dropping  the  negligible  terms,  the  linearized  conservation 
of  momentum  equation  is  found  to  be 

Po^  +  VP  =  0.  (A. 15) 

The  equation  of  state  (A. 8)  is  linearized  by  inspection  of  its  terms,  yielding 


P  =  <%SP- 


(A. 16) 


Using  this  to  eliminate  5p  from  the  linearized  continuity  equation  (A.  13), 

4^r+P°V-u  =  0.  (A. 17) 

c20  dt 

Taking  the  difference  of  the  time  derivative  of  this  equation  with  the  spatial  derivative 
of  the  linearized  momentum  equation  (A.  15)  results  in  the  linear  acoustic  wave  equation 

GS-v2)p=o-  (a-i8) 
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It  is  often  convenient  to  consider  sound  transmission  in  the  frequency  domain  rather 
than  the  time  domain.  The  wave  equation  is  easily  brought  into  the  frequency  domain 
by  applying  a  Fourier  transform,  yielding  the  Helmholtz  equation 


( kl  +  V2)p  =  0, 


(A. 19) 


where  the  acoustic  wavenumber  kQ  =  co/c0  and  u  is  angular  frequency. 

Now  we  will  consider  a  fluid  volume  which  contains  an  acoustic  source.  The  source 
can  be  considered  to  add  a  volume  source  to  the  right  side  of  the  continuity  equation 
(A.  13)  and  a  force  term  to  the  right  side  of  the  momentum  equation  (A.  15).  Carrying 
these  terms  through  the  above  derivation  yields 


(A. 20) 


where  the  volume  and  force  terms  have  been  replaced  by  a  generalized  pressure  source 
JF(x,  t).  This  result  is  known  as  the  inhomogeneous  wave  equation. 

The  radiation  condition  In  general,  solutions  of  the  wave  equation  and  Helmholtz 
equation  can  be  decomposed  into  two  parts  which  respectively  describe  an  outgoing,  or 
forward-traveling  wave  and  an  incoming,  or  backward-traveling  wave.  Let  us  consider  a 
fluid  at  rest.  At  time  t  —  0,  an  acoustic  source  begins  to  generate  sound  in  the  fluid. 
Causality  requires  that  for  t  >  0  the  wave  energy  from  the  source  radiates  outward  and, 
in  the  absence  of  other  sources  or  bodies,  there  is  no  incoming  energy.  The  consideration 
that  all  incoming  solutions  of  the  wave  equation  must  therefore  be  identically  zero  is 
known  as  the  radiation  condition. 
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A. 2  Acoustic  impedance,  intensity,  and  power 

Acoustic  impedance  quantifies  the  sound  absorption  of  a  fluid.  It  is  defined  as  the  ratio 
of  acoustic  pressure  to  the  particle  velocity  in  the  fluid.  For  an  outgoing  wave,  the 
characteristic  impedance  ZQ  of  a  medium  is  given  as 

ZQ  =  -  =  p0c0 ■  (A. 21) 

u 

Rigid  walls  and  pressure  release  surfaces  A  medium  is  considered  acoustically 
hard  or  rigid  if  its  impedance  Z  — >  oo.  A  normally  incident  plane  wave  traveling  through 
a  medium  with  impedance  Z\  is  reflected  back  in  its  entirety  when  it  reaches  a  medium 
of  impedance  Z2  when  Z\  -C  Z2.  If  instead  Z\  Z2l  there  is  again  no  transmission 
and  the  reflected  wave  identically  cancels  the  incident  wave;  the  interface  between  two 
such  media  is  termed  a  pressure  release  surface,  so  named  because  the  acoustic  pressure 
there  is  identically  zero.  Media  whose  impedance  Z  — »  0  are  considered  compliant ,  or 
acoustically  soft. 

Acoustic  intensity  is  the  time  average  of  the  sound  energy  flow  through  a  unit  area 
in  the  direction  normal  to  that  area.  In  integral  form,  acoustic  intensity  I  is  found  as 

1  ptav 

I  =  —  /  pudt.  (A. 22) 

lav  J 0 

The  directivity  of  radiated  sound  is  determined  by  the  dependence  of  intensity  on  the 
direction  in  which  it  is  calculated.  For  sound  which  radiates  omnidirectionally,  intensity 
on  a  spherical  surface  of  radius  r  simplifies,  applying  the  definition  of  impedance,  to  a 
scalar  quantity 

j  =  pO)2 

po^o 

Finally,  the  acoustic  power  W  passing  through  any  surface  S  in  the  fluid  is  the  integral 


(A. 23) 
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of  the  sound  intensity  over  that  surface,  so 


W  =  I  IdS'.  (A. 24) 

Js 

where  S'  is  a  differential  element  of  the  surface  S.  For  omnidirectional  sound,  the  power 
through  a  sphere  of  radius  r  reduces  to 


W  =  Anr2I. 


(A. 25) 
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Appendix  B 

Advanced  topics  from  aeroacoustics 


This  appendix  provides  a  brief  overview  of  select  topics  from  aeroacoustics  pertinent 
to  this  dissertation  and  is  included  for  completeness.  Material  is  taken  from  Howe’s 
Acoustics  of  Fluid- Structure  Interactions  [18],  Theory  of  Vortex  Sound  [20],  and  Hydro¬ 
dynamics  and  Sound  [22]. 

B.l  Acoustic  compactness  and  the  near  and  far  fields 

An  acoustic  source  occupying  a  fluid  region  of  characteristic  dimension  i  and  frequency 
u  radiates  sound  of  wavelength  2nc0/u The  source  is  considered  acoustically  compact 
when  I  is  much  smaller  than  this  wavelength,  or  u)I/c0  =  k0£  <C  1.  The  acoustic  far  field 
exists  at  points  |x|  (taking  the  origin  of  coordinates  at  the  source)  many  wavelengths 
from  the  source,  or  where  fc0|x|  1.  The  near  field  is  thus  that  region  in  the  immediate 
vicinity  of  the  source,  for  which  kQ  <  1. 

Monopoles,  dipoles,  and  quadrupoles  A  monopole  is  an  acoustic  source  which 
can  be  considered  analogous  to  a  pulsating  sphere.  Sound  radiates  from  the  source 
omnidirectionally,  diminishing  in  the  far  field  (as  |x|  — >  cx))  like  A  general  expression 
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for  a  monopole  is  a  source  of  the  form  q(t)5(x),  where  q{t )  is  the  ’strength’  of  the  source. 

A  dipole  is  source  of  the  form  [/j(f)(5(x)]  with  strength  fj(t).  As  a  monopole  can 
be  likened  to  a  pulsating  sphere,  a  dipole  is  analogous  to  an  oscillating  sphere.  Sound 
emanates  from  the  source  as  two  lobes  protruding  in  the  direction  of  that  oscillation  (see 
Figure  B.l),  diminishing  like  jAj  in  the  far  field.  A  point  dipole  is  also  equivalent  to  two 
monopoles  of  opposite  strength  a  short  distance  apart. 


C^O 

dipole 

Figure  B.l.  Directivity  of  the  sound  from  a  dipole  and  quadrupole  source 

A  source  distribution  equivalent  to  four  monopoles  (whose  net  source  strength  is 

d^T' 

zero)  is  called  a  quadrupole.  A  quadrupole  source  has  the  general  form  dx.gx.  (x,  t ),  and 
in  the  far  field  its  generated  sound  diminishes  like  jUp  The  directivity  of  a  quadrupole 
is  depicted  in  Fig.  B.l. 

In  general,  any  acoustic  source  can  be  represented  as  the  sum  of  component  sources  of 
monopole,  dipole,  quadrupole,  and  higher-order  pole  strength.  Typically  it  is  necessary 
only  to  consider  the  leading  order  term  when  determining  the  radiated  acoustic  pressure 
generated  by  the  source. 


quadrupole 
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B.2  Aerodynamic  sound 

The  sound  generated  by  vorticity  in  an  unbounded  fluid  is  called  aerodynamic  sound. 


Lighthill  [32]  developed  the  theory  of  vortex  sound  by  reformulating  the  Navier-Stokes 


equation  into  an  exact,  inhomogeneous  wave  equation  whose  source  terms  are  important 
only  within  the  turbulent  (vortical)  region.  It  was  determined  that  the  solution  in  a 
real  fluid  for  the  density  fluctuations  dp  in  the  far-held  from  a  vortical  acoustic  source 
is  equivalent  to  the  solution  in  an  ideal  fluid  which  is  forced  by  a  stress  distribution  Ttj, 
the  Lighthill  stress  tensor.  This  is  known  as  Lighthill’s  acoustic  analogy.  Including  this 
stress  term  in  the  momentum  equation  used  to  derive  the  acoustic  wave  equation  (see 
Appendix  A)  yields  Lighthill’s  equation-. 


The  problem  of  calculating  aerodynamic  sound  in  a  real  fluid  is  therefore  found  to  be 
equivalent  to  solving  this  equation  for  the  radiation  of  sound  in  an  ideal  fluid  produced 
by  a  distribution  of  quadrupole  sources  whose  strength  per  unit  volume  is  Lighthill’s 
stress  tensor  T^. 


B.3  Acoustic  Green’s  function 

The  Green’s  function  is  an  important  tool  for  solving  inhomogeneous  wave  equations. 


The  free-space  Green’s  function  is  the  outgoing  solution  of  the  wave  equation  generated 


by  the  point  source  5(x  —  y )d(t  —  r)  in  an  unbounded  fluid  located  at  a  point  x  =  y  and 
a  time  t  —  r.  Expressed  mathematically, 
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where  G  =  0  for  t  <  r.  This  has  the  solution 


G(x,  y,  t  -  r)  =  — — ^ - -8  ft-r-  ^  ■  (B.3) 

4vr|x  —  y|  V  co  J 

The  inhomogeneous  wave  equation 

(^-V2)p  =  ^x.<)  (B'4) 

describes  the  behavior  of  waves  generated  by  a  distributed  source  .F(x,  t).  This  source 
can  be  regarded  as  an  array  of  point  sources,  such  that 


JT(y,  r)5(x  -  y )8(t  -  r)d3ydr. 


(B.5) 


The  solution  for  each  constituent  source  of  strength  jF(y,r)h(x  —  y )8(t  —  r)d3ydr  is 
^"(y,  t)G(x,  y,  t  —  r)d3ydr.  Adding  up  these  individual  contributions  we  obtain 


p(x,t) 


JT(y,  t)G(x,  y,  t  -  r)d3ydr 


1 

4vr  _ 

i  r°° 

47 r  /  _ 


-^(y;r) 
»  |x-  y| 

•A  (y,t  - 


8  ( t  —  r 


lx  —  y| 


d3ydr 


lx — y  I 


lx  —  y| 


d3  y. 


(B.6) 

(B.7) 

(B.8) 


Green’s  functions  can  be  applied  to  more  general  acoustic  problems  involving  solid 
or  fluid  boundaries,  and  have  equivalent  expressions  in  both  the  time  and  frequency 
domains.  To  demonstrate  the  use  of  Green’s  functions  in  a  general  problem,  we  will 
consider  a  simple  example  involving  the  solution  of  the  homogeneous  Helmholtz  equation. 
We  consider  an  acoustic  source  on  the  boundary  of  a  volume  of  fluid  defined  by  some 
surface  S.  The  source  generates  an  acoustic  pressure  p(x,uj)  within  the  fluid  volume 
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which  satisfies 

(V2  +  k2a)p  =  0  (B.9) 

subject  to  some  boundary  condition  p  =  p'(x,cv)  on  S.  We  then  introduce  a  Green’s 
function  G(x,u)  with  outgoing  wave  behavior, 

(V2  +  k2)G  =  0  (B.10) 

with  the  boundary  condition  G  =  G'(x,  u)  on  S.  Taking  the  difference  of  G-(B.9)  and 
p-(B.10), 

pV2G-GV2p  =  0.  (B.ll) 

Integrating  over  the  fluid  volume  V, 

[  (pV2G  -  GV2p)  dV  =  0.  (B.12) 

Jv 

Introducing  the  vector  identities 

PV2G  =  V2(pG)  -  V  ■  (GVp)  -  VG  ■  Vp,  (B.13) 

GV2p  =  V2(Gp)  -  V  •  (pVG)  —  Wp  ■  VG  (B.14) 

into  the  integrand  we  find 


[V  •  (pVG)  -  V  ■  (GVp)]  dl/ 


0. 


(B-15) 


Applying  the  divergence  theorem,  this  reduces  to 


[V  ■  (pVG)  -  V  ■  (GVp)]  dV 


dp_ 

dxn 


d^ 


0 


(B.16) 
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where  xn  is  the  vector  outward  normal  to  the  surface  S.  The  acoustic  pressure  p  is  then 
solved  for  by  substitution  into  the  integrand  of  the  Green’s  function  G  and  the  boimdary 
conditions  on  S  specific  to  the  problem’s  geometry  and  source  parameters. 

B.3.1  The  reciprocal  theorem 


XB 


Figure  B.2.  Reciprocal  acoustic  problems  of  equal  point  sources  at  Xa  and  X#. 

An  important  and  very  useful  concept  in  the  application  of  Green’s  functions  for 
solving  acoustic  problems  is  the  reciprocal  theorem.  We  consider  two  acoustic  prob¬ 
lems  in  which  sound  is  generated  in  a  fluid  by  two  point  sources  at  x  =  xa  and 
x  =  xg  in  the  presence  of  a  solid  body  defined  by  a  surface  S.  The  Green’s  func¬ 
tions  G(x,  xa,  uj),  G(x,  xb,co)  which  determine  the  acoustic  pressure  generated  by  each 
respective  source  are  defined  by 

(V2  +  kg)G(x,  xa,  cj)  =  5(x  —  xa)  (B.17) 

(V2  +  fc2)G(x,  xb.w)  =  S(x  —  xB)  (B.18) 
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subject  to  the  boundary  conditions  on  S 


dG  ,  .  G(x,xA,cv)  dG  ,  .  G(x,xB,  uS) 

—  (x,xA,u)  =  -  ,  —  (x,xB,u;)  =  , 

Cs  -Jy  yi  y  .  (jJ  J  (J  Jb  ^  y  5  Wy  J 


(B.19) 


where  xn  is  the  vector  normal  to  the  surface  S  directed  into  the  fluid  and  Z (x,  uj)  is  the 
surface  impedance  on  S.  Mathematically  stated,  the  reciprocal  theorem  holds  that 


G(xa,xb)  =  G(xb,xa), 


(B.20) 


that  is,  the  acoustic  pressure  at  x^  produced  by  the  point  source  at  xB  is  equal  to  the 
acoustic  pressure  at  xB  produced  by  an  equal  point  source  at  X4.  The  reciprocal  theorem 
is  particularly  useful  for  simplifying  acoustic  problems  in  which  the  geometry  local  to 
either  the  source  or  the  point  of  interest  complicates  finding  their  solutions. 


B.3.2  Method  of  images 


y  • 


rigid  wall 

X3  =  0 

y  ° 


Figure  B.3.  Acoustic  image  problem 

Another  important  concept  in  applying  Green’s  functions  is  the  method  of  images. 
Let  us  consider  an  acoustic  source  at  a  point  y  near  an  infinite  rigid  wall  X3  =  0. 
Recognizing  that  the  acoustic  pressure  on  the  wall  vanishes,  we  can  choose  to  reformulate 
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the  problem  by  removing  the  wall  and  introducing  a  source  of  equal  and  opposite  strength 
located  at  the  image  point  y  of  the  original  source  in  the  wall,  that  is,  its  reflection 
through  the  plane  X2  =  0;  the  sound  generated  from  the  image  source  identically  cancels 
that  generated  by  the  original  source  along  X2  =  0,  maintaining  the  boundary  condition 
of  the  original  problem.  This  is  applied  in  the  use  of  Green’s  functions  by  letting  the 
Green’s  function  G  for  the  source  near  the  wall  be  represented  as  the  sum  of  the  free- 
space  Green’s  functions  G'  for  the  sources  at  the  original  point  and  the  image  point, 
or 

G(x,  y,  u)  =  G'(x,  y,  u)  -  G"(x,  y,  u).  (B.21) 
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Appendix  C 


Garfield  Thomas  Water  Tunnel 


Rectangular 
working  section 
508.0  mm  x  114.0  mm 


Figure  C.l.  Schematic  of  the  ARL  0.3048-m  diameter  water  tunnel  facility 

Experiments  using  the  model  snpercavitating  vehicle  are  conducted  inside  ARL’s 
0.3048  m-diameter  water  tunnel  facility.  The  tunnel  is  a  closed  circuit,  closed  jet  system 
which  can  achieve  velocities  up  to  24.38  m/s  under  absolute  pressures  ranging  between 
20.7  and  413.7kPa.  The  flow  is  controlled  by  a  mixed  flow  peerless  pump  driven  by 
a  111.8  kW  motor.  The  tunnel  accomodates  two  test  sections,  a  rectangular  section 
measuring  508mm  by  114.3mm  by  762 mm  long  (shown  above),  and  a  cylindrical  test 
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section  measuring  304.8  mm  in  diameter  by  762  mm  long.  The  cylindrical  test  section  is 
used  for  the  supercavitating  vehicle  tests. 
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Appendix  D 


Equipment  Specifications 


This  appendix  presents  the  pertinent  specifications  for  transducers  and  electronics  used  in 
the  jet  impact  force  experiment.  These  include  the  unsteady  force  sensor,  accelerometer, 
velocimeter,  flowmeter,  signal  conditioner,  and  data  acquisition  hardware. 


D.l  Unsteady  force  transducer 


Model 

Frequency  Response 
Resolution 
Measurement  Range 
Nominal  Sensitivity 


PCB  208C01 
0.01  to  36000  Hz 
0.00045  N-rms 
44.48  N 
112.4  mV /N 


Table  D.l.  Force  transducer  specifications 
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D.2  Accelerometer 


Model 

PCB  352B10 

Frequency  Response 

1  to  17000  Hz 

Resonant  Frequency 

>  65  kHz 

Nominal  Sensitivity 

10  mV/g 

Table  D.2.  Accelerometer  specifications 


D.3  Velocimeter 


Model 

TSI  8345 

Velocity  Range 

0  to  30  m/s 

Temperature  Range 

-17.8  to  93.3°C 

Velocity  Sensitivity 

3%  or  ±0.15  m/s, whichever  is  greater 

Temperature  Sensitivity 

±  0.3°C 

Table  D.3.  VelociCalc  specifications 


This  velocimeter  is  calibrated  annually  by  the  manufacturer. 


D.4  Flowmeter 


Model  Sierra  Instruments  780S-NAA-N5-E2-P3-V1-DD-MP-SR 
Accuracy  ±1% 


Table  D.4.  Flowmeter  specifications 


This  flowmeter  is  calibrated  annually  by  the  manufacturer. 
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D.5  Signal  Conditioner 


Model 

PCB  480B21 

Frequency  Range 

0.15  to  100  kHz 

Voltage  Gain  (±  2%) 

1:1 

(±  2%) 

1:10 

(±  2%) 

1:100 

Spectral  Noise  (1  Hz) (Gain  xl) 

-128  dB 

Spectral  Noise  (10  Hz)  (Gain  xl) 

-140  dB 

Spectral  Noise  (100  Hz) (Gain  xl) 

-144  dB 

Spectral  Noise  (1  kHz)  (Gain  xl) 

-147  dB 

Spectral  Noise  (10  kHz) (Gain  xl) 

-148  dB 

Broadband  Electrical  Noise  (1  to  10  kHz)  (Gain  xl) 

-110  dB/rms 

Spectral  Noise  (1  Hz)  (Gain  xlO) 

-115  dB 

Spectral  Noise  (10  Hz)  (Gain  xlO) 

-118  dB 

Spectral  Noise  (100  Hz)  (Gain  xlO) 

-122  dB 

Spectral  Noise  (1  kHz)  (Gain  xlO) 

-123  dB 

Spectral  Noise  (10  kHz)  (Gain  xlO) 

-129  dB 

Broadband  Electrical  Noise  (1  to  10  kHz)  (Gain  xlO) 

-86  dB/rms 

Table  D.5.  Signal  conditioner  specifications 


D.6  Data  Acquisition  System 


Model  DSP  Technologies  SigLab  50-21 

Frequency  Response  0  to  50  kHz 


Table  D.6.  SigLab  hardware  specifications 


Computer:  Dell  Latitude  D820  laptop  computer 
Software:  Siglab  Dynamic  Signal  Analyzer  (VNA)  v3.2 
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Appendix  E 

Extended  experimental  results  and 
sound  predictions 


This  appendix  presents  all  of  the  unsteady  force  spectra  measured  during  the  experiment 
of  Chapter  4.  Each  test  using  the  parameters  of  injection  rate  Q  and  jet  length  £  (see 
4.1)  was  repeated  three  times;  we  herein  refer  to  these  as  runs  A,  B,  and  C.  Also  included 
are  the  subsequent  predictions  of  the  impingement-generated  sound  power  spectra,  as 
in  §4.5  for  each  data  set.  An  example  of  the  prediction  for  the  self-noise  spectra  at  the 
vehicle  nose  (as  in  §5.5)  is  given  for  each  set  of  jet  parameters  ^  and  Q. 
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E.l  Measured  force  spectra  (§4.4) 


fD0/U0 


Figure  E.l.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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fD0/U0 


Figure  E.2.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  £  =  20,  30,  and  40  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.3.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  £  =  15,  25,  35,  and  45  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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0.001  0.01  0.1  1  10 

fD0/U0 

Figure  E.4.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  £  =  20,  30,  and  40  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.5.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  £  =  15,  25,  35,  and  45  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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Figure  E.6.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  A,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.7.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45  mm,  run  B,  normalized  to  mean  force  Fa  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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fD0/U0 


Figure  E.8.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  £  =  20,  30,  and  40  mm,  run  B,  normalized  to  mean  force  Fa  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.9.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  £  =  15,  25,  35,  and  45  mm,  run  B,  normalized  to  mean  force  F0  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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Figure  E.10.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  B,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.ll.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45  mm,  run  B,  normalized  to  mean  force  Fa  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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fD0/U0 


Figure  E.12.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  B,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.13.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45  mm,  run  C,  normalized  to  mean  force  Fa  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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fD0/U0 


Figure  E.14.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  C,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.15.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45  mm,  run  C,  normalized  to  mean  force  F0  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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Figure  E.16.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0071  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  C,  normalized  to  mean  force  Fa  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 


fD0/U0 


Figure  E.17.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  i  =  15,  25,  35,  and  45nnn,  run  C,  normalized  to  mean  force  Fa  in  15.6  Hz 
bins,  versus  dimensionless  frequency  fD0/U0. 
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fD0/U0 


Figure  E.18.  Measured  force  spectra  for  jets  of  gas  injection  rate  Q  =  0.0094  m3/s  and 
jet  lengths  i  =  20,  30,  and  40  mm,  run  C,  normalized  to  mean  force  F0  in  15.6  Hz  bins, 
versus  dimensionless  frequency  fD0/U0. 
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E.2  Predicted  sound  power  spectra  (§4.5) 


Figure  E.19.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  15  mm,  run  A,  versus  reduced  frequency  kw/D. 
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Figure  E.20.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  20mm,  run  A,  versus  reduced  frequency  kw/D. 


Figure  E.21.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  25  mm,  run  A,  versus  reduced  frequency  kw/D. 
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Figure  E.22.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  30mm,  run  A,  versus  reduced  frequency  kw/D. 


Figure  E.23.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  35  mm,  run  A,  versus  reduced  frequency  kw/D. 
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Figure  E.24.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  40mm,  run  A,  versus  reduced  frequency  kw/D. 
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Figure  E.25.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  45  mm,  run  A,  versus  reduced  frequency  kw/D. 
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Figure  E.26.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  15  mm,  run  B,  versus  reduced  frequency  kw/D. 
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Figure  E.27.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  20  mm,  run  B,  versus  reduced  frequency  kw/D. 
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Figure  E.28.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  25  mm,  run  B,  versus  reduced  frequency  kw/D. 


Figure  E.29.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  30  mm,  run  B,  versus  reduced  frequency  kw/D. 
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Figure  E.30.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  35  mm,  run  B,  versus  reduced  frequency  kw/D. 


Figure  E.31.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  40  mm,  run  B,  versus  reduced  frequency  kw/D. 
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Figure  E.32.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  45  mm,  run  B,  versus  reduced  frequency  kw/D. 


Figure  E.33.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  15  mm,  run  C,  versus  reduced  frequency  kw/D. 
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Figure  E.34.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  20  mm,  run  C,  versus  reduced  frequency  kw/D. 


Figure  E.35.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  25  mm,  run  C,  versus  reduced  frequency  kw/D. 
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Figure  E.36.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  £  =  30mm,  run  C,  versus  reduced  frequency  kw/D. 


Figure  E.37.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  £  =  35  mm,  run  C,  versus  reduced  frequency  kw/D. 
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Figure  E.38.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  40  mm,  run  C,  versus  reduced  frequency  kw/D. 


kw  D 


Figure  E.39.  Predicted  sound  power  spectra  for  jets  of  gas  injection  rate  Q  =  0.0047, 
0.0071,  and  0.0094  m3/s  and  i  =  45  mm,  run  C,  versus  reduced  frequency  kw/D. 
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E.3  Predicted  self-noise  spectra  at  the  vehicle  nose 


(§5.5) 


Figure  E.40.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2 a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
15  mm. 
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Figure  E.41.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
20  mm. 


k  w  a 


Figure  E.42.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
25  mm. 
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Figure  E.43.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
30  mm. 


k  w  a 


Figure  E.44.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
35  mm. 
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Figure  E.45.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
40  mm. 


k  w  a 


Figure  E.46.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0047  m3/s  producing  jets  of  length 
45  mm. 
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Figure  E.47.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
15  mm. 


k  w  a 


Figure  E.48.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
20  mm. 
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Figure  E.49.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
25  mm. 


k  w  a 


Figure  E.50.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
30  mm. 
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Figure  E.51.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
35  mm. 


k  w  a 


Figure  E.52.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
40  mm. 
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Figure  E.53.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0071  m3/s  producing  jets  of  length 
45  mm. 


k  w  a 


Figure  E.54.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
15  mm. 
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Figure  E.55.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
20  mm. 


k  w  a 


Figure  E.56.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
25  mm. 
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Figure  E.57.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
30  mm. 


k  w  a 


Figure  E.58.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
35  mm. 
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Figure  E.59.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
40  mm. 


k  w  a 


Figure  E.60.  Predictions  of  the  sound  pressure  at  the  nose  when  x  =  a,  2a,  3 a  for  an 
ellipsoidal  cavitator  of  aspect  ratio  b/a  =  5/3  for  cases  where  the  cavitator  has  a  pc-nose 
or  is  acoustically  hard  and  gas  is  injected  at  Q  =  0.0094  m3/s  producing  jets  of  length 
45  mm. 
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